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Preface 

The  purpose  of  this  research  twofold.  First,  an 
attempt  was  made  to  validate  the  current  forecasting  method 
utilized  by  HQ  AFLC/DSXR  to  predict  SDT  tonnage.  Secondly, 
if  this  method  was  determined  to  be  statistically  invalid, 
an  effort  was  made  to  develop  a  Box-Jenkins  time  series 
forecasting  model  that  would  produce  more  accurate  and 
reliable  results. 

The  current  method  employed  by  HQ  AFLC/DSXR  was,  in 
fact,  determined  to  be  statistically  invalid  during  the 
course  of  this  research.  Therefore,  an  attempt  was  made  to 
develop  a  Box-Jenkins  model  that  would  offer  more  accurate 
results.  Although  a  significantly  more  accurate  model  was 
not  developed  during  this  research,  the  need  for  developing 
such  a  model  was  given  more  urgency  by  the  fact  that  the 
model  currently  used  is  invalid. 

This  research  would  not  have  been  possible  without  the 
assistance  and  knowledge  of  Mr.  George  T.  Menker  and  his 
staff  at  HQ  AFLC/DSXR.  I  am  also  deeply  indebted  to  my 
advisor  and  friend  Lt  Col  Robert  E.  Trempe  for  his  technical 
knowledge  and  patience  in  assisting  me  with  this  thesis. 
Finally,  I  wish  to  thank  my  wife  Paula  for  her  caring  and 
undying  patience  as  I  struggled  through  this  research  and 
the  entire  AFIT  graduate  program.  Without  her  support  and 
the  brightly  shining  face  of  my  daughter  Allison,  this 
research  would  have  been  much  harder. 


Stephen  L.  Strom 
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Abstract 

This  research  was  conducted  to  analyze  the  Air  Foi ce 
Logistics  Command  (AFLC)  forecasting  method  for  predicting 
Second  Destination  Transportation  (SDT)  tonnage  with  the 
regional  flying  hour  program.  This  thesis  had  two  main 
objectives:  (1)  validate  the  current  forecasting  method  used 
for  computing  tonnage  estimates  to  derive  SDT  budget 
requests,  and  (2)  if  the  current  method's  validity  was  not 
supported,  develop  a  new  forecasting  model,  using  the  same 
input  data,  that  would  produce  more  accurate  and  reliable 
tonnage  estimates. 

Analyzing  graphs  of  the  four  different  data  sets 
researched  in  this  thesis  and  then  conducting  a  statistical 
test  on  the  flying  hour  parameter  for  each  set.  it  was 
determined  that  the  current  method  employed  to  forecast  SDT 
tonnage  was  statistically  invalid  for  two  of  the  four  sets. 
This  determination  was  made  due  to  the  fact  that  the  flying 
hour  parameter  changed  during  the  iterative  regression 
process  used.  This  change  implied  that  SDT  tonnage  and 
flying  hours  were  not  linearly  related. 

Box-Jenkins  (BJ)  time  series  forecasting  models  for 
each  data  set  were  built  and  provided  accurate  and  valid 
forecasts.  For  the  MAC  SDT  time  series,  the  BJ  models  were 
more  accurate  than  the  current  method.  The  BJ  models  for 

vi  i  i 


the  MSC  SDT  time  series,  although  marginally  less  accurate 
than  the  current  method,  were  valid,  whereas  the  current 
method  for  these  two  series  wa3  statistically  invalid. 

This  research  emphasized  the  need  for  an  accurate  model 
and  an  increase  in  the  size  of  the  data  base  used  for 
forecasting.  It  was  also  noted  that  forecasting  SDT  tonnage 
requires  continual  analysis  and  updating  to  ensure  the  model 
being  used  is  appropriate  for  the  data  being  forecasted. 
Finally,  in  statistically  invalidating  the  current  model, 
this  research  has  caused  an  immediate  need  for  an  accurate 
and  valid  model.  Further  research,  particularly  with 
econometric  models,  would  prove  beneficial. 
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AN  ANALYSIS  OF  THE  AIR  FORCE  LOGISTICS  COMMAND  (AFLC) 


FORECASTING  METHOD  FOR  PREDICTING  SECOND 
DESTINATION  TRANSPORTATION  ( SDT ) 

I .  Introduction 

The  Budget  and  Requirements  Division,  Directorate  of 
Programs  and  Resources,  Deputy  Chief  of  Staff  Distribution, 
Headquarters  Air  Force  Logistics  Command  (HQ  AFLC/DSXR),  is 
tasked  with  estimation  and  submittal  of  the  AFLC  budget  for 
Second  Destination  Transportation  (SDT).  These  estimates 
are  derived  from  forecasted  tonnage  requirements  and  are 
submitted  tc  the  Plans  and  Programs  Division,  Directorate  of 
Transportation,  Headquarters  United  States  Air  Force 
(HQ  USAF/LETX ) ,  who,  in  turn,  consolidates  all  SDT  estimates 
throughout  the  Air  Force  into  an  overall  service  SDT  budget 
requirement.  The  major  consequence  of  AFLC/DSXR 
underestimating  its  SDT  requirements  is  that  by  the  end  of 
the  fiscal  year,  SDT  funds  will  be  depleted  and  cargo 
requiring  SDT  will  be  delayed  in  reaching  final 
destinations.  This  is  a  costly  consequence  as  evidenced  by 
a  23  September  1987  letter  from  Brigadier  General  Richard  D. 
Smith,  Deputy  Chief  of  Staff  Materiel  Management, 
Headquarters  Air  Force  Logistics  Command  (HQ  AFLC /MM ) ,  to 
the  Assistant  Secretary  of  the  Air  Force  for  Logistics. 
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The  letter  stated: 

As  requested,  the  cost  of  an  additional  day  of 
shipping  time  of  recoverable  spares  is  $50.8 
million.  With  approximately  24  percent  of  the 
shipments  overseas,  we  estimate  the  cost  of  one 
day  additional  time  at  ports  is  $12.2  million. 

These  numbers  show  the  magnitute  of  impact  of  MAC 
consolidation  of  shipments  at  ports  and  shortfalls 
in  Second  Destination  Transportation  funding 

Another  problem  with  continually  under  or  overesima ting 

is  that  higher  echelons  (i.e.  HQ  USAF/LETX)  lose  confidence 

in  the  method  of  forecasting.  This  latter  problem  has,  at 

times,  plagued  AFLC/DSXR.  Accurate  forecasting  of  SDT 

requirements  will  ensure  all  SDT  shipments  reach 

destination,  and  higher  echelons  will  be  more  apt  to  believe 

in  these  forecasts. 

In  researching  this  topic,  SDT  must  be  clearly  defined. 
To  thoroughly  define  SDT,  the  Air  Force  first  defines  First 
Destination  Transportation  (FDT)  as  the  movement  of  property 
from  the  point  of  origin,  normally  a  contractor  or  supplier, 
to  the  point  where  the  material  is  initially  received  by  the 
Air  Force  for  use  or  storage  for  subsequent  distribution  in 
the  supply  system  (5:124) .  SDT  is  this  subsequent 
distribution  and  includes: 

1.  Port  handling. 

2.  Overocean  transportation. 

3.  Shipments  to  Air  Logistics  Centers  (ALC's)  or 
contractor  depots. 

4.  Shipments  between  bases. 
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5.  Shipments  from  bases  to  repair  facilities  or 
depots.  (5:124;  13;  14:32) 

Currently,  AFLC/DSXR  uses  the  following  linear 
regression  methodology  for  SDT  budget  forecasting. 

Developed  in  the  mid-1970s,  it  proceeds  as  follows: 

1.  The  40  most  recent  quarters  of  historical  tonnage 
and  flying  hour  data  are  obtained,  broken  down  by  overseas 
geographical  region.  The  tonnage  is  also  subdivided  by 
means  of  transport. 

2.  A  series  of  linear  regressions,  beginning  with  the 
initial  40  quarters  of  data,  are  conducted  using  Equation 
(1)  . 


Y  =  0O  +  PiX  (1) 

where:  Y  =  the  dependent,  or  forecasted,  variable,  SDT 

tonnage 

X  =  the  independent  variable,  flying  hours 
Po  =  the  Y-axis  intercept 

(3i  =  the  independent  variable  parameter  value,  also 
the  slope  of  the  regression  line. 


During  each  subsequent  regression,  the  number  of  quarters  of 
data  is  decreased  by  eliminating  the  oldest  quarter  of  data. 
The  final  regression  iteration  uses  the  eight  most  recent 
quarters  . 

3.  For  each  regression,  the  coefficient  of 
correlation  (R) ,  standard  deviation,  and  tonnage  forecasts 
are  computed. 
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4.  The  regression  with  the  highest  R  value  is  chosen 
as  the  equation  to  use,  and  the  resulting  tonnage  forecasts 
are  utilized  to  compute  the  SDT  budget  estimates. 

AFLC/DSXR  also  applies  a  smoothing  technique  to  remove 
the  random  variation  within  the  regression.  After  this 
application,  the  4-step  procedure  above  is  repeated.  The 
final  choice  between  this  latter  technique  and  the  technique 
initially  discussed  above  is  arbitrary.  The  guidelines  for 
directing  this  choice  are  consistency  of  technique  use  over 
a  period  of  time  and  thorough  comparisons  of  the  forecasts 
to  the  most  recent  historical  tonnage  (1:3-4;  12).  In  other 
words,  the  method  selected  for  use  as  the  forecasting  tool 
should  either  consistently  employ  the  nonsmoothing  or 
smoothing  technique.  In  addition,  the  forecasts  should  be 
compared  to  recent  historical  tonnage  to  ensure  there  is  no 
unexplained  change. 

Flying  hours,  the  independent  variable  in  the 
regression,  is  defined  as  the  total  number  of  hours  flown  by 
Air  Force  aircraft  assigned  to  a  specific  overseas 
geographical  region  to  which  each  forecast  applies  as  well 
as  hours  flown  by  transient  aircraft  bedded  down  in  the 
region  in  excess  of  60  days  (12).  These  hours  do  not 
include  those  flown  by  the  MAC  aircraft  hauling  SDT  to  these 
regions . 

The  forecasted  requirements  are  divided  into  five 
means  of  transport: 
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1.  Military  Airlift  Command  (MAC). 

2.  Military  Sealift  Command  (MSC). 

3.  Military  Traffic  Management  Command  (MTMC),  which 
includes  port  handling  operations. 

4.  Logistics  airlift  (LOGAIR). 

5.  Government  Bills  of  Lading  (GBL's),  which  include 
commercial  air  and  surface  transportation  (13;  14:32,34). 
These  are  the  five  predominant  means  of  transport  utilized 
in  shipping  SDT  within  the  Air  Force. 

The  Air  Force  flying  hour  program  by  overseas 
geographical  area  is  used  as  the  independent  variable  in 
forecasting  future  tonnage  requirements,  the  dependent 
variable,  for  each  means  of  transportation  except  LOGAIR 
(10:3;  12;  14:34).  This  functional  relation  is  shown  in 
Equation  (1)  above.  These  means  of  transport,  with  the 
exception  of  LOGAIR,  support  five  overseas  geographical 
areas  (major  commands):  European  (USAFE) ,  Pacific  (PACAF), 
Alaskan  (AAC),  Southern  (USAFSO) ,  and  Northeastern  (northern 
radar  sites)  (10:3;  12). 

The  forecasted  tonnage  and  subsequent  SDT  budget 
estimates  are  computed  separately  by  geographical  region  for 
the  following  reasons: 

1.  Costs  of  cargo  movement  to  and  within  each  region 
differ  due  to  the  weight  of  the  shipment  and  the  distance 
moved . 
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2.  Flying  hours  for  aircraft  assigned  to  each  region 
are  projected  by  the  major  command  (MAJCOM)  associated  with 
that  region  (10:4).  Hence,  several  different  combinations 
of  SDT  tonnage  forecasts  are  computed. 

Table  1  depicts  AFLC  tonnage  and  cost  for  each  mean  in 
FY  88. 


Table  1 


Tonnage  and  Cost  by  Means 


Fiscal  Year  1988  (13) 


Means 


Tonnage  (000) 


Cost  (millions) 


MAC  83,563  S/T*  $130 
MSC  1,098,820  M/T**  68 
MTMC  984,774  M/T  14 
LOGAIR  65,675  S/T  76 
GBL  189,545  S/T  43 
TOTAL  $331 


*  A  short  ton  (S/T)  is  equal  to  2,000  pounds. 


**  A  measurement  ton  (M/T)  is  equivalent  to 
approximately  40  cubic  feet.  (20:838) 


As  noted  in  the  Table  1 ,  MAC  and  MSC  shipments  were 
allotted  a  total  of  $198  million,  or  57  percent  of  AFLC 1 s 
SDT  funding. 

The  SDT  cost  by  theater  is  also  an  important 
consideration  in  the  budgeting  process.  This  cost  breakdown 
allows  AFLC  to  determine  where  the  majority  of  its  SDT 
expenditures  occur.  Table  2  lists  the  percent  of  SDT 
funding  by  MAJCOM  based  on  FY  89  figures. 
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Table  2 


Percent  of  AFLC  SDT  Cost  in  FY  89  by  MAJCOM  (13) 

MA J COM  Percent  of  Cost 

USAFE  37 

CONUS  29 

PACAF  25 

Northern  3 

USAFSO  3 

AAC  2 

Other  1 

As  shown  in  Table  2,  SDT  shipments  to  USAFE  and  PACAF 
account  for  62  percent  of  AFLC 1 s  FY  89  SDT  funding. 

Shipments  to  these  same  two  theaters  account  for  87  percent 
of  AFLC ' s  SDT  funding  for  overseas  areas. 

Specific  Probl cm 

An  accurate  and  reliable  forecasting  method  to  predict 
SDT  requirements  is  necessary  to  ensure  proper  budget 
submittals.  A  reliable  forecasting  method  is  free  of  random 
error;  therefore,  when  forecasting  with  a  reliable  method, 
the  forecasts  will  show  little  or  no  random  error,  or 
differences  from  the  actual  values  (6:98).  An  accurate 
forecast  refers  to  the  ability  of  the  method  to  predict 
future  values.  Accuracy  is  a  comparative  measurement  and 
the  method  exhibiting  the  best  accuracy  measurement  is 
usually  selected  for  forecasting  (11:567).  Several  methods 
for  measuring  forecast  accuracy  are  discussed  in  Chapter  II. 
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The  current  forecasting  method  utilized  by  AFLC/DSXR 
for  computing  tonnage  estimates  to  derive  future  SDT  budget 
requests  was  developed  in  the  mid-1970s  and  has  not  been 
validated  in  recent  years.  Validation  of  this  current 
method  would  support  or  fail  to  support  its  accuracy  and 
reliability  as  a  forecasting  tool.  In  validating  the  model, 
evidence  must  be  found  and  presented  showing  that  the  method 
forecasts  what  it  is  suppose  to  forecast.  In  addition, 
evidence  must  be  presented  supporting  a  causal  relationship 
between  the  variables  in  the  forecasting  model  (4:37;  6:94). 

Justification 

The  SDT  budget  estimates  derived  from  the  forecasted 
tonnage  requirements  are  submitted  to  the  Plans  and  Programs 
Division,  Directorate  of  Transportation,  HQ  USAF/LETX,  who, 
in  turn,  consolidates  all  SDT  estimates  throughout  the  Air 
Force  into  a  service  SDT  budget  requirement  (12).  Recently, 
USAF/LETX  expressed  concern  with  respect  to  the 
effectiveness  of  the  current  forecasting  method  in 
predicting  future  tonnage  requirements  (18).  This  concern 
stems  from  the  cutbacks  in  SDT  as  well  as  other 
transportation  funding.  As  a  result  of  the  recent  funding 
cuts  throughout  the  armed  services,  it  is  more  critical  than 
ever  to  be  able  to  predict,  with  a  greater  amount  of 
accuracy,  future  tonnage  and  funding  requirements.  Accurate 
predictions  are  necessary  to  ensure  SDT  funding  is  available 
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throughout  the  fiscal  year.  Table  3  lays  out  the  actual  SDT 


funding  for  AFLC  over  the  past  five  fiscal  years. 


Table  3 

AFLC  SDT  Funding  (millions)  (13) 


FY 

Requirement 

Fundinq 

Difference 

85 

$407 

$407 

$0 

86 

407 

407 

0 

87 

410 

410 

0 

88 

385 

315 

-70 

89 

460 

337 

-73 

As  noted  in 

Table  3,  AFLC 1 s 

SDT  Program 

has  been  funded 

less  than  100  percent  over 

the  past  two 

fiscal  years. 

Scope 

Since  87  percent  of  overseas  SDT  occurs  within  USAFE 
and  PACAF  by  MAC  and  MSC  movements,  AFLC/DSXR  believes  that 
if  the  forecasts  for  these  areas  and  transportation  means 
are  validated,  its  forecasting  method  is  reliable  and 
accurate  (12).  For  this  reason,  the  research  for  this 
thesis  was  limited  to  the  examination  of  the  forecasts  for 
these  two  areas  by  these  two  means  of  transport.  In 
addition,  this  research  was  limited  to  the  use  of  flying 
hours  by  overseas  geographical  location  (USAFE  or  PACAF)  as 
the  independent  variable  to  predict  the  dependent  variable, 
SDT  tonnage.  This  latter  limitation  was  made  to  offer 
AFLC/DSXR  the  best  possible  forecasting  method  for  the  data 


it  currently  uses,  and  also  due  to  a  lack  of  da*  ■?.  for  other 
possible  independent  variables. 

Research  Objectives 

This  research  was  conducted  to  assist  in  validating  the 
forecasting  method  utilized  by  AFLC/DSXR  to  predict  SDT 
tonnage.  These  forecasts  must  be  accurate  to  ensure  proper 
budget  requests  are  submitted  for  approval.  Therefore,  the 
research  objectives  of  this  thesis  were  to: 

1.  Validate  the  current  forecasting  method  used  for 
computing  tonnage  estimates  to  derive  future  SDT  budget 
requests . 

2.  If  the  current  method’s  validity  was  not 
supported,  develop  a  new  forecasting  model,  using  the  same 
input  data,  that  would  produce  more  accurate  and  reliable 
tonnage  estimates. 

Resear ch  Hypothesis 

The  research,  or  null,  hypothesis  (Ho)  for  this  thesis 
is  that  the  forecasting  method  utilized  by  AFLC/DSXR  to 
compute  SDT  tonnage  estimates  is  an  accurate  and  reliable 
predictor  of  future  SDT  budget  requirements.  The  alternate 
hypothesis  (H^)  is  that  the  method  used  is  not  an  accurate 
and  reliable  predictor. 

In  research,  only  two  conditions  (or  states  of  nature) 
exist,  these  being  the  null  hypothesis  is  true  or  it  is 
false.  In  addition,  only  one  of  two  possible  decisions 
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about  these  states  of  nature  is  made  --  the  null  hypothesis 
is  either  accepted  as  true  or  it  is  rejected  as  false.  The 
combination  of  these  two  conditions  and  two  decisions  leads 
to  four  possible  situations.  Two  situations  lead  to  correct 
decisions,  while  the  remaining  two  result  in  erroneous 
decisions  (6:352-353).  See  Figure  1. 

Condition 

Hq  -  true  Hq  -  false 

Decision 

Accept  Hq  Correct  Decision  Type  II  error 

Confidence  level  probability  =  3 

probability  =  1-a 

Reject  Hq  Type  I  error 

Significance  level 
probability  =  a 

Figure  1.  Hypothesis  Testing 
Condition  -  Decision  rlatrix  (6:353) 

Type  I  error  (a)  occurs  when  a  true  hypothesis  is 
rejected  and  Type  II  error  ((3)  occurs  when  a  false 
hypothesis  is  accepted  (6:353). 

In  research,  the  problem  with  accepting  the  null 
hypothesis  is  that  it  cannot  be  proven,  but  only  supported 
or  not  supported.  One  reason  for  this  problem  is  that  there 
always  remains  a  possibility  that  statistics  will  not  detect 
a  difference  between  the  null  and  alternate  hypotheses  when, 
in  fact,  there  was  a  difference.  Secondly,  future  results, 
using  more  powerful  statistical  tests  or  analyses  with 


Correct  Decision 
probability  =  1-3 
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higher  correlation  between  variables,  cannot  be  known  in 
advance  (4:44) . 

Although  the  null  hypothesis  cannot  be  proved  or 
disproved,  decisions  must  be  made  as  though  the  null 
hypothesis  is  true  or  false.  Acceptance  or  rejection  of  the 
null  hypothesis  should  therefore  be  based  on  the  confidence 
estimated  in  relation  to  established  standards  (4:45). 

Several  factors  must  be  considered  when  setting 
standards  for  this  research.  First  the  cost  of  implementing 
a  new  computer  forecasting  software  package  is  of  major 
concern.  Along  with  a  new  system  comes  training,  and  this 
factor  is  not  costless  and  must  therefore  be  considered. 
Increased  confidence  in  the  forecasts  by  higher  echelons  is 
also  a  factor.  The  new  method  should  be  considered  for 
possible  implementation  if  it  provides  more  accurate  and 
reliable  forecasts  and  higher  echelons  are  more  confident  in 
the  results . 

For  this  research,  the  confidence  level  probability,  or 
the  probability  of  accepting  the  null  hypothesis  when  it  is 
true,  was  set  at  95  percent.  In  other  words,  the 
probability  of  making  a  Type  I  error  (a)  was  set  at  5 
percent.  This  probability  and  associated  null  hypothesis  is 
considered  "conservative."  This  conservativeness  implies 
that  if  the  null  hypothesis  is  rejected  as  false,  there  is 
almost  no  chance  that  the  null  hypothesis  was  in  fact  true. 
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AFLC/DSXR,  in  reviewing  this  research,  should  select 
confidence  level  and  Type  I  error  probabilities  that  meet 
its  concerns  based  on  the  factors  it  considers  relevant. 
However,  the  confidence  level  and  Type  I  probability  error 
selected  for  this  research  provide  adequate  protection  to 
AFLC/DSXR  fjom  incurring  unnecessary  switching  costs  while 
subjecting  the  current  forecasting  methodology  to  a  rigorous 
test  of  its  accuracy. 

Plan  of  Analysi s 

The  remaining  chapters  of  this  thesis  are  devoted  to 
the  actual  development  and  implementation  of  the  research. 
Chapter  II  is  devoted  to  a  review  of  the  literature 
concerning  the  current  forecasting  method  as  well  as  the 
general  topic  of  forecasting.  Chapter  III  lays  out  the 
methodology  used  to  conduc t  this  research  in  pursuit  of  the 
objectives  stated  earlier.  The  results  and  analysis  of  this 
research  are  presented  in  Chapter  IV.  Finally,  Chapter  V  is 
devoted  to  the  findings  and  implications  of  this  research. 
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Literature  Review 


II  . 

This  literature  review  is  composed  of  three  parts.  The 
first  section  briefly  explains  the  concept  of  forecasting 
and  offers  an  outline  of  forecasting  methods  as  well  as  a 
group  of  methods  for  measuring  forecast  accuracy.  The 
second  part  details  the  current  forecasting  method  utilized 
by  AFLC/DSXR  in  forecasting  future  SDT  tonnage.  Finally,  in 
the  third  section,  previous  reports  and  studies  conducted  on 
AFLC/DSXR' s  current  forecasting  method  are  reviewed  and 
analyzed . 

Fo  r ecasting 

Forecasting  is  conducted  for  two  main  purposes: 
extrapolation  and  intervention  analysis.  Extrapolation  is 
predicting  how  future  events  will  develop  based  on  past- 
events.  Intervention  analysis  predicts  the  effect  of 
management  decisions  or  policy  changes,  and  if  an  effectiv-- 
change  does  occur,  assists  in  determining  the  direction  and 
magnitude  of  the  change  (2:12.18-12.19). 

Forecasting  was  conducted  in  this  thesis  for  the 
purpose  of  extrapolation.  Predictions  of  future  SDT  tonnage- 
are  needed  to  ensure  proper  budget  submittals.  Caution  T,usf 
be  exercised  with  regard  L  O  CA  trapolating  past  the  scope  of 
the  data  used  to  construct  the  model.  However,  in  order  r 
provide  a  forecast  of  the  future,  extrapolation  must  be 
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conducted.  Accuracy  in  forecasting,  therefore,  is 


essential,  and  without  it,  resultant  decisions  based  on 
these  forecasts  will  be  erroneous. 

Typology  of  Forecasting.  There  are  many  different 
categorizations,  or  typologies,  of  forecasting  methods; 
however,  most  typologies  include  the  same  methods  and  only 
the  groupings  differ.  There  are  five  main  categories  of 
forecasting  methods  including: 

1.  Subjective. 

2 .  Ad  hoc . 

3.  Causal  or  structural. 

4.  Deterministic. 

5 .  Time  series  . 

Subjective .  Subjective  forecasting  methods  are 
heavily  reliant  on  the  opinions  of  the  forecaster  and 
involve  techniques  such  as  an  educated  guess,  engineering 
judgment,  or  Delphi.  Subjective  models  are  intuitive  and 
are  the  most  commonly  used  forecasting  technique.  Reasons 
for  using  a  subjective  approach  are  short  response  times, 
costs  of  using  a  more  sophisticated  technique,  and  lack  of 
historical  data.  The  main  characteristic  of  subjective 
forecasts  is  that  they  cannot  be  replicated  due  to  the  lack 
of  an  actual  "model"  (2:12.7-12.8;  3). 

Ad  Hoc.  Ad  hoc  forecasting  methods  are  based 
primarily  on  past  history  and  work  well  in  areas  of  apparent 
stability.  Ad  hoc  models  include,  but  are  not  limited  to: 
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1.  Simple  average. 

2.  Simple  moving  average. 

3.  Weighted  moving  average. 

4.  Exponential  smoothing. 

5.  Holt’s  Two-Parameter  Model. 

6.  Winter's  Three-Parameter  Model. 

The  simple  average  is  accurate  when  there  is  stability, 
or  stationari ty ,  within  the  data  accompanied  by  little 
variation.  Stationarity  implies  there  is  no  trend  in  the 
data  or  that  the  mean  of  the  data  remains  constant.  The 
simple  moving  average,  or  moving  average,  eliminates  the 
need  to  store  large  amounts  of  data  necessary  to  use  the 
simple  average;  however,  valuable  information  maybe  lost  if 
old,  but  informative,  data  are  excluded.  The  weighted 
moving  average  allows  the  forecast  to  give  differing  weights 
to  different  observations.  This  alleviates  the  problem  in 
using  the  moving  average  where  all  observations  are  given 
equal  weights.  Exponential  smoothing  is  useful  if  the  data 
are  stationary.  Holt's  Two-Parameter  Model  works  well  on 
data  with  both  stationary  and  trend  components.  A  trend 
occurs  when  the  variance  of  the  data  remains  the  same,  but 
the  mean  changes  up  or  down.  Finally,  if  the  data  exhibit  a 
seasonal  component  as  well  as  stationary  and  trend 
components,  the  Winter's  Three-Parameter  Model  works  best. 
Seasonality  is  defined  as  periodicity  within  the  data  and 
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may  occur  every  12  periods  within  monthly  data  or  every  four 
periods  within  quarterly  data  (2:12.8-12.10,12.16;  3;  11:14- 
16)  . 

Causal .  Causal,  or  structural,  forecasting 
methods  are  often  referred  to  as  econometric  models.  The 
forecast  is  a  set  of  functions  based  on  assumed  or 
theoretical  causal  relationships  in  the  data.  Simple  linear 
and  multiple  regression  models  are  of  this  type.  An  example 
of  a  simple  linear  regression  model  is  in  Equation  12). 
Simple  linear  regression  is  a  special  case  of  multiple 
regression  in  that  only  one  independent  variable  is  used. 


Y  =  0Q  +  0iX  +  e  (2) 

where:  Y  =  a  series  of  historical  values  of  a  dependent 

variable 

X  =  a  corresponding,  or  causal ,  series  values  of  an 
independent  variable 
Po  =  the  Y-axis  intercept 

3l  =  independent  variable  parameter  value;  slope  of 
the  line 
e  =  error  term. 


A  multiple  regression  model  example  is  in  Equation  (3) . 


Yi  -  Po  +  3lxli  +  3 2X2  i  +  -  -  •  +  Pkxki  +  e  <3) 

where:  Y ±  =  dependent  variable  observation 

xki  =  the  observed  value  of  the  independent  variable 
Pk  =  the  independent  variable  parameter  value 
e  =  error  term. 


AFLC/DSXR's  current  method  is  similar  to  equation  (2) 
with  Y  equal  to  SDT  tonnage  and  X  equal  to  flying  hours.  In 
causal  model  development,  the  forecaster  must  ensure  all 


pertinent  variables  are  included  and  nonexplanatory 
variables  are  excluded  (2:12.12-12.14;  3;  9:10-11). 

Two  problems  involved  with  regression  are 
heteroscedas tici ty  and  autocorrelation .  Heteroscedasticity , 
or  unequal  variances  within  the  data,  does  not  normally 
occur  in  time  series  studies.  This  is  due  to  the  likelihood 
of  the  changes  in  the  dependent  variable  and  the  changes  in 
any  number  of  the  independent  variables  being  of  the  same 
magnitude.  In  regression,  however,  variances  of  the  error 
term  are  not  used  to  prove  the  absence  of  bias.  Therefore, 
the  problem  encountered  in  regression  when 

heteroscedasticity  is  present  is  that  biased  estimates  of 
parameter  variances  are  used  thus  the  statistical  tests  and 
confidence  intervals  used  are  incorrect  (17:141-142). 

Autocorrelation,  or  serial  correlation,  occurs  when 
errors  associated  with  a  given  observation  carry  over  into 
other  periods.  Generally,  autocorrelation  does  not  affect 
the  unbiasedness  of  regression  model  estimators,  but  does 
affect  the  estimators'  efficiency.  Under  positive 
autocorrelation,  the  parameter  estimates  are  concluded  to  be 
more  precise  than  they  actually  are,  and  this  leads  to  the 
tendency  to  reject  the  null  hypothesis  when,  in  fact,  it  is 
true  (17 : 152-153) . 

Deterministic.  Deterministic  forecasting  methods 
make  the  dependent,  or  observed,  variable  a  function  of  time 
and  are  of  the  form  in  Equation  (4) . 
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f  (t) 


(4) 


Yt  = 

where:  f(t)  =  a  function  of  time. 

This  type  of  model  forecasts  the  future  behavior  of  a  time 
series  based  on  past  behavior.  Deterministic  implies  that 
reference  is  not  made  to  the  sources  or  causes  of  randomness 
within  the  data  (17:473).  As  in  causal  models, 
heteroscedasticity  and  autocorrelation  may  present  problems 
and  affect  the  ultimate  accuracy  of  the  model. 

Time  Series.  Time  series  models  identify  patterns 
in  the  history  of  a  variable  and  use  the  information 
gathered  to  predict  future  values  of  that  same  variable. 
These  models  are  used  when  there  is  no  knowledge  about  the 
causal  relationships  affecting  the  variable  being  forecasted 
(17:470-471).  Most  time  series  methods  assume  the  following 
characteristics: 

1.  The  time  period  at  which  each  data  value  is 
measured  is  of  equal  length  (i.e.  years,  quarters,  months). 

2.  The  data  values  are  arranged  in  time  sequence  from 
the  earliest  to  latest  time  period,  and  there  are  no  missing 
values . 

3.  The  process  and  method  of  measurement  are 
consistent  over  time  (9:5). 

Time  series  forecasting  models  follow  one  of  two 
approaches:  self-pro jecting  or  cause-and-ef feet  (9:9). 


19 


Self-projecting,  or  univariate,  methods  are  based  on 
the  history  of  the  times  series,  and  range  in  sophistication 
from  moving  averages  to  Box-Jenkins.  These  methods  are 
useful  for  several  reasons: 

1.  They  are  quick  and  easy,  and  can  handle  large 
numbers  of  series  in  an  efficient  manner. 

2.  They  require  only  a  small  amount  of  data  and  are 
normally  inexpensive  to  perform. 

3.  They  provide  accurate  short-  to  mid-term  forecasts. 
Self -pro jecting  methods  do  not,  however,  account  for  any 
external  factors  that  may  influence  the  time  series  being 
forecasted  (9:9-10). 

Cause-and-ef f ect  methods  of  time  series  forecasting 
account  for  external  influences  by  the  use  of  mathematical 
relationships  between  the  series  being  forecasted  and  any 
number  of  other  series  representing  the  influencing  factors. 
These  methods  are  useful  because: 

1.  More  information  is  available  concerning  the 
forecast . 

2.  Interrelationships  between  many  factors  are  taken 
into  account . 

3.  Accurate  forecasts  can  be  provided  for  the  mid-  to 
long-term . 

Cause-and-ef feet  methods  of  time  series  forecasting  require 
larger  amounts  of  data  and  usually  require  more  time  and 
money  to  perform  compared  to  other  methods  (9:10-11). 
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The  time  series  model  used  in  this  thesis  and  explained 
further  below  is  the  Box-Jenkins  model.  Box-Jenkins  (BJ) 
models  use  past  values  of  the  data  and/or  error  terms  in  the 
model  development.  A  simplistic  example  of  a  BJ  model  is  in 
Equation  ( 5 ) . 


Yt  =  Ft  +  Et  (5) 

where:  =  time  series 

Ft  =  any  patterns 

Et  =  random  error  term.  (2:12.16;  9:17) 

The  BJ  method  involves  four  phases  in  the  model 
building  process,  including: 

1.  Pattern  identification. 

2.  Model  specification. 

3.  Model  diagnosis. 

4.  Hypothesis  testing  and  forecasting. 

In  building  a  forecasting  model  for  a  particular  time 
series,  the  process  above  becomes  an  iterative  one.  The 
final  model  selection  from  those  specified  in  this  iterative 
process  is  based  on  comparisons  of  the  diagnostic  checks 
conducted  in  the  third  step.  The  following  text  more 
clearly  defines  each  of  these  four  steps. 

First,  the  underlying  patterns  in  the  raw  data  must  be 
identified  (i.e.  trend,  seasonality,  cycle).  Visual 
inspection  of  the  data  plot  allows  possible  identification 
of  an  array  of  patterns.  A  pattern,  opposite  of  randomness, 
is  defined  as  a  reoccurrence  of  data.  One  type  of  pattern 
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is  known  as  trend.  A  trend  occurs  when  the  variance  of  the 
data  remains  the  same,  but  the  mean  may  change  either  up  or 
down.  If  the  mean  of  the  data  remains  constant,  then  no 
trend  exists  and  the  data  are  referred  to  as  stationary. 
Seasonality,  a  second  pattern,  is  defined  as  periodicity 
within  the  data.  Seasonality  may  occur  every  12  periods 
within  monthly  data  or  every  four  periods  within  quarterly 
data.  Another  pattern  is  cyclical  in  nature  and  is 
identified  by  a  consistently  rising  or  falling  pattern  over 
time.  Cyclical  and  seasonal  patterns  differ  in  that  a 
cyclical  pattern  normally  includes  two  or  more  seasonal 
ones.  (2:12.16;  3;  9:46;  11:14-16). 

Next,  the  appropriate  model  must  be  specified  based  on 
the  identified  patterns.  In  specifying  the  appropriate  BJ 
model  based  on  the  patterns  identified,  there  are  an  array 
of  choices.  The  basi''  RT  mode),  assuming  stationarity ,  has 
two  basic  parameters:  autoregressive  parameters  (AR)  and 
moving  average  parameters  (MA)  (9:48-49). 

The  AR  parameters  use  past  values,  or  patterns,  of  the 
variable  to  predict  future  values  as  shown  in  Equation  (6). 

Y  t  =  $lYt_i  +...+  $pYt_p  +  Et  (6) 

where:  Yp  =  the  time  series 

$p  =  the  AR  parameter 
p  =  order  of  the  model 
Yp-p  =  past  values  of  the  time  series 
Ep  =  the  error  term.  (9:49-50) 
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The  HA  parameters  use  the  random  errors,  or  error 
terms,  occurring  in  past  time  periods  to  predict  the  time 
series  value  at  time  t  as  shown  in  Equation  (7) .  MA 
forecasting  models  do  not  depend  on  past  values  of  the  time 
series  in  computing  forecasts. 

Yj-  =  -©3_E]_  — .  .  QqE^-q  +  E^  (7) 

where:  Y<-  =  the  time  series 

0q  =  the  MA  parameter 

q  ~  order  of  the  model 

Et-q  ~  the  error  term.  (9:51) 

When  specifying  a  model  based  on  the  patterns 
identified,  AR  and  MA  parameters  may  be  included  in  the  same 
model.  These  types  of  models  are  referred  to  as  ARMA  models 
and  follow  the  format  in  Equation  (8). 

Yt  =  *lYt-l  *pYt_p  -  (G1Et_1  +  ...h  0qEt_q)  +  Et  (3) 

where:  Y^  =  the  time  series 

3>p  =  the  AR  parameters 

Yt-p  =  past  values  of  the  time  series 
0q  =  the  MA  parameters 

Et-q  =  past  values  of  the  error  term.  (9:52) 

If  the  identified  patterns  reveal  the  time  series  to  be 
nonstationary,  differencing  of  the  data  eliminates  this 
nonstationarity .  Differencing  a  series  involves  subtracting 
the  first  value  of  the  series  from  the  second,  the  second 
from  the  third,  and  so  on  (11:36).  ARMA  models  that  include 
differencing  are  referred  to  as  autoregressive,  integrated, 
moving  average  ( ARIMA )  models  (Equation  (9)). 
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Zt  =  ♦iZt  — i  >...■*■  ^pZ^-p  _  (02.Et-i  +.  .  .+  0qE^_q)  +  E^  (9) 

where:  Zt  =  the  differenced  (stationary)  time  series 

<t»p  =  the  AR  parameters 

Zp-p  =  past  values  of  the  differenced  time  series 
0q  =  the  MA  parameters 
Et-q  =  past  values  of  the  error  term 

Et  =  the  error  term.  (9:129,138;  17:529-538) 

The  third  step  is  to  diagnose  the  model  using  certain 
"tools"  or  computer  outputs.  These  tools  are  explained  in 
the  methodology. 

Finally,  after  the  model  passes  the  diagnostic  checks, 
the  forecast  must  be  conducted  (2:12.16-12.17;  3;  9:47; 
11:251-254) . 

Methods  of  Measuring  Forecast  Accuracy .  Several 
methods  exist  to  measure  the  accuracy  of  forecasting 
methods,  and  the  resultant  measurement  values  can  be  used  to 
compare  different  forecasting  methods.  The  accuracy 
measurement  method  used  is  important  because  it  may 
determine  which  forecasting  method  is  actually  used  to 
forecast  into  the  future  (2:12.17).  Several  accuracy 
measurement  methods  exist  including: 

1.  Average  error  (AE) . 

2.  Mean  absolute  deviation  (MAD) . 

3.  Mean  square  error  (MSE). 

4.  Standard  deviation  (SD). 

5.  Mean  absolute  percent  error  (MAPE) .  (19:113-115) 
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Many  of  these  methods  are  based  on  the  error  of  the 
forecast.  This  error  is  many  times  referred  to  as  the  error 
term,  forecast  error  or  bias,  random  error  component,  or 
randomness.  No  matter  what  it  is  referred  to  as,  the 
random  error,  or  error  term,  is  almost  always  denoted  by  E 
or  e  in  any  equation  or  formula.  The  error  term  is  defined 
as : 

The  part  of  a  series  that  has  no  explanation 
outside  of  chance  alone.  That  is,  any  data 
collected  to  represent  some  process  or  activity 
will  always  be  in  error  owing  to  measurement 
errors  or  outside  random  influences.  Each  data 
value  in  a  time  series,  therefore,  is  expected  to 
be  in  error  from  the  "true"  value  for  that  period. 

The  manner  in  which  these  errors  occur  is 
generally  assumed  to  be  random.  (If  they  weren’t, 
they  would  demonstrate  some  relationship  and 
therefore  be  part  of  the  pattern  component.)  (9:271) 

In  regression,  as  well  as  any  other  forecasting  model 

employing  the  error  term,  basic  assumptions  about  the  term 

are  usually  made.  These  assumptions  are: 

1.  The  has  a  normal  distribution  with  a  mean  of  zero 
and  a  standard  deviation  of  o. 

2.  The  error  terms  for  each  time  period  are 
independent  of  one  another  (16:339). 

The  average  error  (AE)  is  the  sum  of  the  forecast 
errors  divided  by  the  number  of  forecast  periods.  The  AE , 
expressed  mathematically  in  Equation  (10),  usually 
understates  the  magnitude  of  the  forecast  error;  however, 
the  sign  of  the  AE  helps  to  reveal  the  direction  of  forecast 
error  (19:113) . 
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AE  =  (E]_  +  E2  +  .  .  .  +  En)  /  n  (10) 

where:  En  =  the  error  in  the  n*-*1  time  period 
n  =  number  of  periods  in  the  series. 

The  mean  absolute  deviation  (MAD)  provides  an  accurate 
measurement  of  the  magnitude  of  the  forecast  error.  The 
MAD,  expressed  in  Equation  (11),  is  the  sum  of  the  absolute 
values  of  the  forecast  errors  divided  by  the  number  of 
forecast  periods. 

MAD  =(! EX !+  1 E2 !+...+! En !)/  n  (11) 

where:  !  En )  =  absolute  value  of  the  error  at  time  period  n 

n  =  number  of  periods  in  the  series. 

While  the  MAD  is  an  accurate  measure  of  the  forecast  error, 
the  use  of  absolute  values  in  the  MAD  removes  any  evidence 
of  forecast  bias  (19:113). 

As  expressed  in  Equation  (12) ,  the  mean  square  error 
(MSE)  is  the  sum  of  the  squared  forecast  error  values 
divided  by  the  number  of  forecast  periods  (19:113). 

MSE  =  (Ex2  +  E22  +  .  .  .  +  En  2 )  /  n  (12) 

where:  En2  =  the  squared  value  of  the  error  at  period  n 
n  =  number  of  periods  in  the  series. 

The  MSE  gives  more  importance  to  large  errors  in  its 
computation  than  small  errors.  The  MSE  also  removes  any 
evidence  of  forecast  bias  (11:19). 
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The  standard  deviation  (SD) ,  shown  in  Equation  (13),  is 
the  square  root  of  the  sum  of  the  squared  forecast  error 
values  divided  by  the  number  of  forecast  periods  minus  one 
( 19 : 115) . 

SD  =  [(EX*  +  E22  +  .  .  .  +  En2)  /  (n  -  1)]1/2  (13) 

where:  En 2  =  the  squared  value  of  the  error  at  period  n 
n  =  number  of  periods  in  the  series. 

The  SD  does  not  give  an  indication  as  to  the  direction  of 
the  forecast  error. 

Finally,  the  mean  absolute  percent  error  (MAPE)  finds 
the  absolute  percentage  error  for  each  individual  forecast 
and  computes  the  average  of  these  forecasts  as  a  whole.  The 
MAPE,  shown  in  Equation  (14),  has  the  advantage  of  allowing 
comparisons  between  different  time  series  which  is  not 
permitted  or  possible  when  using  the  MSE  (11:19)  .  This 
advantage  is  due  to  the  fact  that  the  MAPE  expresses  the 
forecast  error  as  a  percentage  of  the  mean  of  the  data  being 
forecasted . 

MAPE  =  f  c  (  ;  EX  :  /Y!  >  +  .  .  .  +  OEni/Yn)]  /  nl  100%  ill' 

where:  En  =  the  error  at  period  n 

n  =  number  of  periods  in  the  series 
Yn  =  the  actual  tonnage  at  period  n. 

In  using  any  of  these  accuracy  measurement  methods  to 
compare  a  group  of  forecast  models,  the  model  exhibiting  the 
smallest  accuracy  measurement  value  can  be  said  to  be  the 
most  accurate  in  the  group. 
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Current  AFLC/DSXR  Forecasti nq  Method ology 

As  stated  previously  in  Chapter  I,  AFLC/DSXR  uses  a 
linear  regression  methodology  for  SDT  budget  forecasting. 
Developed  in  the  mid-1970s,  it  proceeds  as  follows: 

1.  The  40  most  recent  quarters  of  historical  tonnage 
and  flying  hour  data  are  obtained,  broken  down  by  overseas 
geographical  region.  The  tonnage  is  also  subdivided  by 
means  of  transport. 

2.  A  series  of  linear  regressions,  beginning  with  the 
initial  40  quarters  of  data,  are  conducted  using  Equation 
(15)  . 

Y  =  P0  +  Plx  '  15  ' 

where:  Y  =  the  dependent,  or  forecasted,  variable  SDT 

tonnage 

X  =  the  independent  variable  flying  hours 
(3 q  =  the  Y-axis  intercept 

P>1  =  the  independent  variable  parameter  value,  or  the 
slope  of  the  line. 

During  each  subsequent  regression,  the  number  of  quarters  of 
data  is  decreased  by  eliminating  the  oldest  quarter  of  data. 
The  final  regression  iteration  uses  the  eight  most  recent 
quarters . 

3.  For  each  regression,  the  coefficient  of 
correlation  (R) ,  standard  deviation,  and  tonnage  forecasts 
are  computed. 

4.  The  regression  with  the  highest  R  value  is  chosen 
as  the  equation  to  use,  and  the  resulting  tonnage  forecasts 
are  utilized  to  compute  the  SDT  budget  estimates. 
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AFLC/DSXR  also  applies  a  smoothing  technique  to  remove 
the  random  variation  within  the  data.  After  this 
application,  the  4-step  procedure  above  is  repeated.  The 
final  choice  between  this  latter  technique  and  the  technique 
initially  discussed  above  is  arbitrary.  The  guidelines  for 
directing  this  choice  are  consistency  of  technique  use  over 
a  period  of  time  and  thorough  comparisons  of  the  forecasts 
to  the  most  recent  historical  tonnage  (1:3-4;  12) . 

Previous  Reports  and  Studies 

The  need  for  an  accurate  predictor  of  tonnage  to  use  in 
computing  budget  requirements  is  not  new.  In  fact,  over  the 
past  several  years,  many  reports  have  been  written 
concerning  this  subject  (1;  7;  8;  10).  The  reports  and 
studies  reviewed  in  this  chapter  were  conducted  on  the 
method  currently  utilized  by  AFLC/DSXR  to  predict  SDT 
tonnage . 

Foster  Study  (1977).  During  an  SDT  review  conducted  by 
the  Office  of  the  Secretary  of  Defense  (OSD)  in  1976,  the 
reviewer  raised  concern  about  the  validity  of  AFLC ' s 
forecasting  method  for  predicting  SDT  tonnage  and  its 
reliarce  on  the  flying  hour  program  (7:1).  This  study, 
conducted  and  written  by  Newton  W.  Foster  for  the 
Directorate  of  Management  Sciences,  Deputy  Chief  of  Staff 
Plans  and  Programs,  Headquarters  Air  Force  Logistics  Command 
(HQ  AFLC/XRS) ,  was  in  direct  response  to  the  OSD  review. 
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The  objectives  of  the  study  were: 

1.  To  support  the  use  of  flying  hours  as  a  predictor 
of  SDT . 

2.  To  develop  a  better  method  of  predicting  SDT  if  the 
flying  hour  related  computation  could  not  be  supported 
(7:1)  . 

Regression  analysis  was  used  in  conducting  this  study 
and  the  data  used  as  variables  in  this  technique  were 
relatable  to  transportation.  Sixteen  quarters  of  data,  from 
FY  73/1  through  FY  76/4,  were  collected  and  categorized 
under  21  different  headings  (i.e.  manpower,  requisitions, 
overseas  flying  hours,  worldwide  flying  hours).  These  data 
were  collected  by  six  major  geographical  regions:  PACAF , 
USAFE,  AAC ,  USAFSO,  Northeastern,  and  Worldwide.  The 
means  of  transport  within  each  region  in  which  data  were 
collected  were:  MAC,  MSC ,  and  GBL  (worldwide  only)  (7:3-4). 

Two  types  of  regression  were  used:  simple  linear 
regression  and  multiple  regression.  Functionally,  these 
regressions  are  expressed  in  Equation  (16). 

Y  =  f (Xi)  ( 16 ) 

where:  Y  =  SDT  tonnage,  the  dependent  variable 

Xj_  =  the  independent  variable;  for  simple  linear 

regression  -  any  one  of  the  21  categories,  and 
for  multiple  regression  -  any  combination  of  the 
21  categories  . 

These  analyses  were  conducted  on  each  mean  of 
transport  for  each  geographical  region  using  the  21 
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categories  of  data  as  independent  variables.  Simple  linear 
regression  involved  two  variables:  SDT,  the  dependent 
variable;  and  each  of  the  21  data  categories  taken  one  at  a 
time  as  the  independent  variable.  Multiple  regression 
involved  the  use  of  more  than  one  independent  variable  and 
the  dependent  variable  (7:4-5). 

The  conclusions  of  this  study  were: 

1.  The  forecast  method,  although  not  a  totally  valid 
predictor  of  SDT  tonnage,  was  the  most  logical  predictor 
based  on  the  data  provided  and  examined. 

2.  A  better  forecasting  method  for  predicting  SDT 
tonnage  for  a  particular  geographical  region  by  a  specific 
transportation  mode  was  not  evident  based  on  the  data 
provided  and  examined  (7:6). 

Grayson  Research  Study  (1977).  Major  John  Grayson,  an 

Air  Command  and  Staff  College  (ACSC)  student,  conducted 

research  for  partial  fulfillment  of  the  requirements  for 

ACSC  graduation.  Conducting  this  research  as  a  result  of 

its  direct  utility  to  the  Air  Force,  Grayson  stated: 

Increasingly  austere  funding  has  challenged  the 
Air  Force  to  maintain,  and  even  improve,  operational 
readiness  of  the  fleet  with  fewer  available  spares. 

One  step  toward  meeting  this  challenge  has  been 
increased  reliance  on  the  transportation  network 
to  provide  more  timely  distribution,  or  redistribution 
of  available  spares  from  supply  to  user  locations. 

(8:1) 

He  further  noted  that  funding  for  this  transportation 
network,  including  SDT,  had  long  been  a  problem  within  the 
Department  of  Defense.  The  purpose  of  this  study  was  to 


investigate  AFLC ' s  method  for  determining  SDT  requirements 
and  ultimately  offer  recommendations  to  improve  the  process 
(8:4)  . 

In  reviewing  the  AFLC  forecasting  method  and  its 
reliance  on  the  Air  Force  flying  hour  program  (FHP) ,  Grayson 
concluded  that  "...there  is  an  overreliance  on  the  FHP  in 
determining  SDT  requirements"  (8:24-25).  Grayson  theorized 
that  using  flying  hours  as  the  sole  independent  variable  to 
predict  SDT  (shown  in  Equation  (17))  assumed  a  statistical 
insignificance  of  many  other  variables  impacting  SDT. 

Y  =  f(X>  (17) 

where:  Y  =  SDT  tonnage,  the  dependent  variable 

X  =  regional  flying  hours,  the  independent  variable. 

Two  examples  of  other  potentially  significant  variables  were 
manpower  and  differing  amounts  of  logistical  support  per 
flying  hour  for  different  types  of  aircraft  (8:17,19). 

Grayson,  as  a  result  of  this  research,  recommended  two 
actions  be  taken.  First,  efforts  should  continue  by  the  Air 
Staff  and  AFLC  to  identify  programs  whose  cargo  can  be 
separated  from  the  general  cargo  classif ication .  Since  the 
SDT  tonnage  included  in  the  general  cargo  category  is 
predicted  from  the  FHP,  the  separation  of  other  cargo 
classifications  would  reduce  the  overall  percentage  of  SDT 
tonnage  dependent  on  the  FHP.  Secondly,  it  was  recommended 
that  AFLC  identify  other  potential  factors  influencing  SDT 
and  develop  a  multivariate  model  (8:25-26). 
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Lamb-Sarnacki  Thesis  (1978).  Captains  Christopher  J. 
Lamb  and  Joseph  B.  Sarnacki,  Air  Force  Institute  of 
Technology  (AFIT)  graduate  students,  stated  the  problem  was 
the  lack  of  a  valid  forecasting  method  for  predicting  SDT 
tonnage  and,  thus,  an  inability  to  compile  accurate  SDT 
budget  estimates.  The  objectives  of  this  thesis  were  to: 

1.  Determine  the  significance  of  flying  hours  and 
manpower  as  reliable  predictors  of  SDT. 

2.  Develop  a  computer  model  to  predict  future  SDT 
tonnage  if  either  or  both  of  these  variables  were  proved 
reliable  (10:1,10). 

The  forecasting  technique  used  in  this  thesis  was 
discontinuous  linear  regression  ( DLR ) .  This  technique  is 
used  whenever  a  linear  regression  function  changes  slope  and 
at  the  same  time  shifts  vertically  up  or  down,  and  the 
researchers  believed  this  to  be  the  case.  An  additional 
variable  was  needed  to  account  for  this  vertical  shift  in 
the  data  (10:15;  15:316).  Equation  (18),  quoted  directly 
from  the  thesis,  functionally  defines  the  general  form  of 
the  DLR  model  used. 

Z  =  Bq  +  BiX  +  fb2(X-X1)  Xpl  +  B3Xd1I* 

+  ( B  4 ( X~X2 ) Xp2  +  B5XD2| * 

+  (B6(X-X3)  Xp3  +  B7XD3I* 

+  B8Y  +  f  B g ( Y-Yx / Ypi  +  B10Ydi I  * 

+  lB11(Y-Y2)Yp2  +  B12YD2  I* 

+  (B13(Y-Y3)Yp3  +  B ! 4 Yd3 I  *  +  e  (18) 
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*  Discontinuity  adjustments 

where:  Z  =  Tonnage  transported  by  MAC, 

X  =  Flying  hours  (either  programmed  or  actual), 

Y  =  Manpower  (either  programmed  or  actual), 

Xi ,  X2 ,  X3,  Y 1 ,  Y2,  Y3  =  Points  where  data  becomes 

discontinuous , 

Xpi ,  Xp2 ,  Xp3,  Ypi,  Yp2 ,  Yp3, 

XD1'  XD2'  xD3 '  yD1'  yD2'  yD3  =  Indicator  (dummy 

variables)  defined  as: 

xpl  =  XD1  =  I'  if  x  >  X1 >  otherwise  0 

xp2  =  XD2  =  I-  if  X  >  X2;  otherwise  0 

xp3  =  XD3  =  1.  if  X  >  X3;  otherwise  0 

Ypi  =  Ypi  =1,  if  Y  >  Yj;  otherwise  0 

Yp2  =  yD2  =  1 '  if  Y  >  Y2;  otherwise  0 

Yp3  =  yD3  =  1>  if  Y  >  Y3;  otherwise  0 

Bg,...,Bi4=  Coefficients  of  regression,  and 

e  =  The  explained  component  in  each  value 

that  is  not  explained  by  the  independent 
variables.  (10:18-19) 

Only  MAC  data  were  evaluated  and  separate  models  were 
derived  for  each  of  five  geographical  regions  (European, 
Pacific,  Alaskan,  Southern,  and  Northeastern)  as  well  as  the 
total  database.  The  thesis  concluded  that  the  models 
developed  were  more  accurate  than  the  current  method 
employed  by  AFLC  based  on  the  mean  absolute  deviation  (MAD) 
(10:38) . 
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Several  advantages  were  noted  for  using  the  method 
developed  in  this  thesis  including: 

1.  The  lack  of  necessity  for  an  arbitrary  base  period 
for  each  area  since  SDT  estimates  would  now  relate  to  direct 
chanjec  \ r  f  1  v’  ricj  ho’iirs  ?.nd  !”3rr 

2.  Faster  results  could  be  obtained  from  the  use  of  a 
computer-assisted  software  package  such  as  the  Statistical 
Package  for  the  Social  Sciences  (SPSS)  used  in  this  1973 
AFIT  thesis. 

3.  A  more  accurate  and  valid  estimate  of  quarterly  SDT 
tonnage  would  be  provided. 

Abell  Study  (1982).  In  1982,  Joseph  A.  Abell  conducted 

an  evaluation  of  AFLC ' s  forecasting  method.  Abell  stated: 

The  rationale  for  the  use  of  flying  hours  as  an 
indicator  of  tonnage  movements  is  a  result  of  past 
experience.  That  is,  it  has  been  shown  that  in 
most  of  the  five  geographical  areas  served  by  MAC 
and  MSC ,  the  variation  in  the  trends  in  tonnage 
movements  can  be  related  to  the  variation  in  the 
number  of  hours  flown  in  the  specific  geographical 
area.  If  it  is  assumed  that  past  relationships 
will  remain  fairly  constant,  then  it  can  also  be 
assumed  that  one  variable,  programmed  flying 
hours,  can  be  used  to  predict  another,  tons.  (1:1) 

Initially,  Abell  believed  duplication  of  prior 

forecasts  was  necessary  to  properly  evaluate  AFLC ' s 

forecasting  method;  however,  as  a  result  of  the  inability  to 

duplicate  the  forecasts,  a  study  evolved  to  determine  the 

reason  for  this  disagreement.  An  evaluation  of  the  data 

smoothing  technique  employed  by  AFLC  was  also  conducted  to 

determine  its  effect  on  the  data  and  forecast  accuracy  (1:5) . 
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Abell  used  the  CREATE  system  and  the  Statistical 
Package  for  the  Social  Sciences  (SPSS)  to: 

1.  Develop  scattergrams  of  the  data  used  in  the 
regressions . 

2.  Estimate  the  regression  equations  and  accompanying 
statistics. 

3.  Compute  necessary  statistics  on  each  of  the 
individual  variables  (1:5). 

Simple  linear  regressions  were  conducted  on  all  five 
areas  (Europe,  Pacific,  Alaska,  Southern,  and  Northeastern). 
An  evaluation  of  the  October  1981  forecasts  was  also 
conducted  (1:6) . 

Throughout  his  research,  Abell  was  unable  to  duplicate 
the  results  of  the  regressions  being  evaluated.  However,  he 
stated: 

The  nonduplication  of  prior  results,  again,  does 
not  mean  that  the  use  of  linear  regression  is 
invalid.  To  me  it  meant  that  there  was  a  "kink" 
in  computer  program  currently  used  and  that  it 
could  be  corrected.  (1:6) 

Abell  used  three  checks  for  computational  accuracy 
applicable  to  linear  regression.  First,  the  regression  line 
must  pass  through  a  point  (X,Y),  where  X  is  equal  to  the  sum 
of  the  flying  hours  for  a  particular  region  divided  by  the 
total  number  of  observations,  or: 

X  =  (Xi  +  X2  +  ...  +  Xi)  /  n  (19) 

where:  Xn  =  flying  hours  at  period  i  (where  i  =  l,2,...n) 

n  =  the  total  number  of  observations. 
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Secondly,  by  substituting  X  into  the  developed  regression 
equation,  Y  becomes:  Y  =  a  +  bX  .  Any  deviation  from  the 
actual  value  of  Y  would  then  be  considered  forecast  error. 
The  final  check  involved  substituting  the  actual  X  value 
into  the  developed  regression  equation  and  subtracting  this 
result  from  the  actual  Y  value.  The  sum  of  these 
differences  should  be  zero.  Abell  determined  AFLC 
forecasting  method  was  devoid  of  computational  accuracy. 
This  discovery  led  to  search  for  the  cause  (1:6-8) . 

Abell  assumed  as  each  new  regression  was  conducted,  it 
was  done  so  independent  of  all  other  regressions.  He  also 
assumed  the  data  used  for  each  regression  would  produce  the 
same  results  --  or  simply  stated,  the  results  were 
replicatable .  Abell,  however,  found  these  assumptions  fals 
(1:8). 

Abell  concluded  there  was  a  "kink"  in  AFLC ' s  program 
that  added  together  the  total  number  of  observations  from 
all  previous  regressions  with  the  number  of  observations  in 
the  next  to  last  regression.  For  example,  if  the  first 
regression  was  run  with  40  observations,  the  second 
regression  used  40  +  39  =  79  observations;  40  +  39  +  38 
117  observations  were  used  in  the  third,  and  so  on.  The 
number  of  observations  should  have  decreased  by  one  during 
each  successive  regression.  For  example,  if  40  observations 
were  used  in  the  initial  regression,  then  39  would  be  used 
in  the  second,  38  in  the  third,  until  the  last  regression 
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was  conducted  with  the  eight  most  recent  observations. 

Abell  believed  this  "kink"  generated  forecasts  incapable  of 

being  duplicated  (1:9). 

Abell  offered  the  following  conclusion: 

The  use  of  regression  analysis  as  a  method  to 
forecast  future  tc"-r._gc  movements  over  MSC  end  Mar 
should  be  continued  and  is  appropriate  for  the 
following  reasons: 

1)  it  is  dependable  and  defendable 

2)  it  xs  able  to  incorporate  the  effects  of 
past  trends  into  the  estimate 

3)  it  is  able  to  incorporate  indicators  of 
increases/decreases  of  future  operations 
into  the  estimate 

4)  provides  a  measure  of  the  probable  error 
in  the  estimate 

5)  provides  a  measure  of  the  strength  of  the 
relationship  between  tonnage  movements 
and  flying  hours,  the  correlation 
coefficient.  (1:22-23) 

Interpretation  of  Previous  Reports  and  Studies.  It  can 
be  inferred  from  these  reports  and  evaluations  that  there  is 
disagreement  as  to  the  validity  of  AFLC/DSXR's  forecasting 
method.  The  reports  reviewed  in  this  chapter  have  shown  the 
current  method  utilized  is  suspect  and  therefore 
questionable  as  a  forecasting  tool.  However,  only  the 
Lamb-Sarnacki  thesis  offered  a  solid  alternative  to  the 
current  method,  but  it  was  never  implemented.  One 
observation,  in  reviewing  these  sources,  is  that  no  attempt 
has  been  made  to  use  more  sophisticated  forecasting 
techniques  such  as  ad  hoc  or  time  series  methods.  These 
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techniques  include  Holt's  Two  Parameter,  Winter' 
Parameter,  and  Box-Jenkins  methods. 


Three 
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Ill . 


Methodology 


The  primary  objective  of  this  research  was  to  validate 
the  current  forecasting  method  used  by  AFLC/DSXR  to  compute 
SDT  tonnage  estimates.  If  the  method's  validity  as  a 
forecasting  tool  could  not  be  supported,  an  attempt  would  be 
made  to  develop  a  valid  model  using  the  same  input  data. 

This  chapter  explains  the  methodology  followed  in  conducting 
this  research. 

Data  Acquisition 

To  initiate  this  research,  the  necessary  flying  hour 
and  SDT  tonnage  data  were  obtained  from  AFLC/DSXR.  This 
office  gathers  this  data  in  order  to  accomplish  its  current 
forecasting.  Historical  quarteily  data  for  flying  hours  and 
SDT  tonnage  were  obtained  for  FY  78/1  through  FY  37/4.  As 
explained  in  Chapter  I  of  this  thesis,  shipments  by  MAC  and 
MSC  to  PACAF  and  USAFE  account  for  87  percent  of  AFLC ’ s  SDT 
funding  for  overseas  areas.  Actual  quarterly  SDT  tonnage 
data  were  also  collected  for  FY  88/1  through  FY  89/2  in 
order  to  compare  the  forecasts  for  these  same  quarters. 

These  data  are  included  in  Appendix  A. 

Current  Model  Val idation 

In  conducting  this  research,  no  exact  method  could  be 
found  to  evaluate  the  type  of  model  AFLC/DSXR  utilizes; 
however,  after  reading  literature  on  this  subject,  it  was 
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determined  that  if  the  flying  hour  parameter  (|3*  in  equation 
(20))  in  the  model  could  be  shown  to  change  statistically 
after  any  iteration  was  conducted,  the  model  was 
structurally  unstable.  In  using  this  iterative  linear 
regression  approach,  AFLC/DSXR  assumes  that  p]_  does  not 
change  and,  therefore,  that  the  data  are  stationary. 

However,  if  does  change,  this  iterative  approach  becomes 
invalid . 

Y  =  (30  +  BiX  (20) 

where:  Y  =  the  dependent,  or  forecasted,  variable  SDT 

tonnage 

X  =  the  independent  variable  flying  hours 
|3q  =  the  Y-axis  intercept 

=  the  flying  hour  parameter,  or  the  slope  of 
the  line. 

If  this  instability  was  present,  it  implied  the  model  was 
invalid.  In  addition,  this  instability  would  indicate 
nonstationarity  within  the  data,  and  therefore  the  iterative 
linear  regression  approach  would  not  be  useful  in  predicting 
future  values  of  the  time  series. 

In  an  attempt  to  determine  the  stability  of  the  flying 
hour  parameter  as  successive  iterations  of  the  forecast 
mode]  were  conducted,  the  following  steps  were  conducted. 

1.  The  flying  hour  parameter,  P x ,  for  each  iteration 
was  computed.  These  values  were  computed  using  QUATTRO,  a 
spreadsheet  software  package  developed  by  Borland 
International . 
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2.  The  standard  error  (sg^)  for  each  flying  hoar 
parameter  was  computed  with  the  use  of  QUATTRO. 

3.  Using  these  standard  errors,  95  percent  confidence 

intervals  for  ($]_  were  established  for  each  iteration  using 
equation  (21)  (15:58-60;  16:492-493). 

Pl±,:^/2s[n  (  i1 ) 

where:  ta/2  =  the  value  of  the  two-tailed  test-static-ic 

for  a  =  .05  and  n  -  2  degrees  of  freedom. 

4.  If  any  two  of  these  confidence  intervals  did  not 
overlap,  then  the  model  was  showed  structurally  unstable, 
and  was,  therefore,  invalid  as  a  forecasting  tool. 

The  null  hypothesis  (Hq)  for  validating  the  current 
model  was  that  the  flying  hour  parameters  (3i  in  equation 
(20))  for  each  regression  iteration  were  equal.  Hq  thus 
implied  the  current  forecasting  model  used  by  AFLC/'DSXR  was 
stable  and  valid. 

The  alternate  hypothesis  (H^)  was  that  at  least  one  of 
the  flying  hour  parameters  of  any  regression  iteration  was 
not  equal  to  the  others.  therefore  implied  the  model  was 

structurally  unstable  and  therefore  invalid.  If  this  was 
the  case,  it  was  also  inferred  that  the  data  was 
nonstationary . 

The  decision  rule,  or  test  statistic,  shown  in 
equation  (21) ,  was  used  to  draw  inferences  about  the 
parameters.  Hq  was  rejected  if  the  decision  rule  proved 
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that  any  one  of  the  flying  hour  parameters  was  not  equal  t 
the  others;  otherwise,  the  Hq  was  not  rejected. 


The  hypothesis  testing  was  conducted  based  on  the 
following  format; 


H0 :  Pi , 8  =  Pi , 9  =  • • •  =  Pi, 40 

HA:  3l,i  /  Pi, 3  =  =  P 1 , n 

where:  i  =  any  one  iteration  conducted  with  3 
to  40  periods  of  data 
n  =  total  number  of  iterations 
conducted  excluding  i 

Test  Statistic:  P^  ±  ta/2Spi 

where:  -a/2  =  the  value  of  the  test-statistic  f 

a  =  .05  and  n  -  2  degrees  of 
freedom  (df ) . 

Rejection  Region:  Reject  Hq  if  any  two  of  the 

regression  iteration  confidence 
intervals  did  not  overlap. 


QUATTRO  was  used  for  conducting  this  hypothesis 
testing . 

Mod_e  1_  App  1  i  c  a  t  i  o  n 

SDT  tonnage  forecasts  were  obtained  using  AFLC/DSXR's 
forecasting  method,  where  a  series  of  linear  regressions, 
beginning  with  the  last  40  quarters  ( FY  73/1  through  o7/4) 
of  historical  data,  were  conducted  using  equation  (22; . 


Y  =  P0  +  PxX 


where : 


Y 

X 


Po 

Pi 


the  historical  Sht  tonnage 
the  historical  flying  hours 
the  y-axis  intercept 

the  flying  hour  parameter,  or  the  slope. 


During  each  subsequent  regression,  the  number  of 
quarters  of  historical  data  was  decreased  by  eliminating  the 
oldest  quarter  of  data.  The  last  regression  was  conducted 
utilizing  the  eight  most  recent  quarters  of  historical  data. 
The  model  producing  the  highest  coefficient  of  correlation 
(R)  was  selected  as  the  forecasting  model.  R  is  described 
functionally  in  equation  (23). 

R  =  ±  [1  -  ( SSE  /  SSTO)]1/2  (23) 

where:  SSE  =  the  variation  in  the  dependent  variable  (Y) , 

or  SDT  tonnage,  when  the  independent 
variable  (X),  or  flying  hours,  is  used  in 
the  regression  model 

SSTO  =  the  variation  in  Y  when  X  is  not  taken  into 
account.  (15:89-90) 

Model  Building 

Predicated  on  instability,  Box-Jenkins  (BJ)  time  series 
modeling  was  chose. i  as  the  method  to  use  in  developing  a  new 
model  to  forecast  SDT  tonnage.  This  choice  was  made  because 
BJ  identifies  patterns  in  the  history  of  a  time  series  and 
uses  these  patterns  to  develop  the  appropriate  model. 

In  conducting  the  BJ  model  building  and  forecasting 
process  a  computer  software  package  was  used.  TIMES,  mainly 
written  at  The  Ohio  State  University,  was  used  for  this  time 
series  analysis  and  forecasting.  The  necessary  cutput  to 
evaluate  and  select  the  final  models  chosen  was  also 
obtained  using  TIMES. 

In  building  a  Box-Jenkins  time  series  model,  four 
distinct  steps  were  followed. 
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1.  A  plot  of  the  raw  data  as  well  as  other  statistics 
(i.e.  the  autocorrelation  function,  ACF,  and  the  partial 
autocorrelation  function,  PACF)  were  compiled  and  reviewed 
to  identify  any  underlying  patterns  within  the  data  (i.e. 
trend,  seasonality,  etc.).  The  ACF  uses  statistical 
measures  to  determine  the  strength  of  correlation  of  time 
series  values  at  a  certain  number  of  periods  apart.  It 
determines  the  relation  of  a  time  series  to  itself  over 
time.  The  PACF  also  uses  statistical  measures  to  determine 
the  strength  of  other  relations  within  the  time  series 
(9:54-55).  These  ACFs  and  PACFs  were  compiled  using  TIMES. 
In  evaluating  these  functions,  the  values  at  lags  one 
through  three  were  considered  significant  if  they  were 
greater  than  one-half  of  the  corresponding  standard  error. 
The  values  at  lags  four  through  six  were  considered 
significant  if  they  were  larger  than  the  standard  error.  If 
the  values  at  lags  seven  through  nine  were  greater  than  two 
times  the  corresponding  standard  error,  they  were  considered 
significant.  Finally,  all  the  remaining  lagged  values  were 
significant  if  they  were  at  least  three  times  as  large  as 
the  corresponding  standard  error. 

2.  Based  on  these  patterns,  an  appropriate  ARMA  or 
ARIMA  model  was  specified  with  estimated  AR  and  MA 
parameters  from  the  40  quarters  of  time  series  data. 

3.  Diagnostic  tests  were  conducted  on  the  model  to 
ensure  it  was  an  accurate  and  appropriate  model.  There  were 
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eight  tests  conducted  including: 

a.  The  residuals  were  plotted  and  reviewed. 
Ideally,  the  mean  of  the  residuals  should  be  zero  or  near 
zero.  There  should  also  be  no  change  in  the  variance  of  the 
residuals.  If  these  conditions  are  not  met,  there  is  a 
problem  in  the  model  in  that  all  the  patterns  within  the 
time  series  have  not  been  accounted  for  by  the  specified 
model . 

b.  The  residual  autocorrelation  function  (ACF) 
was  inspected.  If  there  are  any  significant  values,  or 
"spikes",  at  any  of  the  lags,  the  model  does  not  explain  all 
the  variation  in  the  data. 

c.  The  Portmanteau  Lack  of  Fit  Test  value,  or  the 
Q-statistic,  was  evaluated.  Ideally,  the  value  should  be 
small  compared  to  other  Q-statistics  with  the  same  number  of 
degrees  of  freedom.  If  it  becomes  too  large,  the  model  does 
riot  fit  the  data  (9:93-94;  17:549-550). 

d.  The  cumulative  periodogram  of  the  estimated 
residuals,  or  errors,  was  evaluated.  The  periodogram  should 
be  linear  or  very  close  to  linear,  otherwise,  the  residual 
variation  in  the  model  cannot  be  attributed  to  purely  random 
processes  . 

e.  A  histogram  of  the  residuals  was  reviewed. 

The  residuals  should  be  normally  distributed  about  zero  with 
a  constant  variance.  If  the  histogram  was  mound-shaped  and 
centered  about  zero,  this  assumption  was  considered  true. 
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f.  The  sum  of  the  squared  errors  (SSE)  and  the 
mean  squared  errors  (MSE)  of  each  time  series  were 
evaluated.  Ideally,  these  values  should  be  small.  Large 
values  indicate  the  model  was  not  a  good  predictor  of  the 
time  series. 

g.  The  Fourier  Transform  of  the  autocorrelations, 
or  Power  Spectrum,  was  evaluated.  This  graoh  should  be 
smooth  and  horizontal,  and  anything  to  the  contrary  implies 
a  problem  in  the  model. 

h.  The  Schwarz  statistic,  or  the  Bayesian 
Information  Criterion  (BIC),  was  reviewed.  Equation  (24) 
was  used  to  compute  t'hp  BT^  which  assesses  goodness  of  fit 
as  well  as  penalizing  for  complexity  within  the  model.  This 
complexity  is  concerned  with  the  number  of  parameters  in  the 
model.  The  BIC,  ideally  small  in  size,  normally  leads  to 
models  of  lower  orders  than  other  similar  statistical  tests. 

BIC  =  SSE  *  n(p+q)/n  (24) 

where:  SSE  =  the  sum  of  the  squared  errors 

n  =  the  number  of  periods  in  the  data 

p  =  the  number  of  autoregressive  parameters 

q  =  the  number  of  moving  average  parameters. 

As  seen  from  Equation  (24),  if  two  models  have  nearly  the 
same  SSE,  but  one  of  the  models  has  more  parameters  than  the 
other,  the  model  with  fewer  parameters  would  be  chosen. 

4.  Finally,  after  a  model  was  built  that  passed  these 
tests  with  stability,  the  forecast  was  conducted  (3;  9:90- 
97;  17:543-552).  The  forecasts  were  conducted  using  TIMES 
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and  the  models  developed  during  this  building  process. 

An  important  point  here  is  that  this  4-step  process  is 
iterative  in  nature  and  therefore  more  than  one  model  is 
usually  specified  before  a  final  choice  is  made.  This 
allows  for  a  comparison  of  the  diagnostics  between  the 
models.  For  this  research,  only  the  final  models  for  each 
of  the  four  data  sets  will  be  discussed. 

Forecast  Accuracy  Measurement 

As  explained  in  Chapter  II,  there  are  several  accuracy 
measurement  methods  that  can  be  used  to  compare  the  accuracy 
of  forecasts  between  different  models.  The  choice  of 
accuracy  measurement  method  is  important  because  it  may 
determine  which  forecasting  method  is  ultimately  used. 

The  mean  absolute  percent  error  (MAPE) ,  shown  in 
equation  (25) ,  was  used  to  compare  the  forecasts  of 
AFLC/DSXR ' s  method  and  the  models  built  during  this 
research . 

MAPE  =  f  [  (  ,*  EX !  /YX )  +  .  .  .  +  (|En:/Yn)]  /  nl  100%  (25) 

where:  En  =  the  error  at  period  n 

n  =  number  of  periods  in  the  series 
Yn  =  the  actual  tonnage  at  period  n. 

The  forecast  errors  for  AFLC/DSXR's  models  were 
computed  by  subtracting  the  actual  SDT  tonnage  from  the 
forecasts.  Likewise,  the  errors  for  the  BJ  models  developed 
in  this  research  were  computed  by  subtracting  the  actual  SDT 
tonnage  from  the  predictions. 
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In  comparing  AFLC/DSXR's  model  and  the  models  developed  in 
this  thesis,  the  forecast  model  producing  the  smaller  MAPE 
was  deemed  to  be  more  accurate. 

Summary 

This  chapter  laid  out  the  steps  used  to  conduct  this 
research.  The  research  began  with  data  acquisition  followed 
by  a  statistical  test  to  validate,  or  invalidate,  the 
current  model  as  a  forecasting  tool.  Next,  the  forecasts  of 
the  current  model  were  verified  by  quickly  calculating  these 
values  using  the  AFLC/DSXR  method.  Following  these 
calculations,  BJ  models  were  built  to  forecast  the  four  data 
sets  being  researched.  Finally,  the  AFLC/DSXR  models  were 
compared  to  the  BJ  models  for  forecast  accuracy.  This 
methodology  leads  into  the  results  and  analysis  discussed  in 


Chapter  IV. 
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This  chapter  includes  results  obtained  during  research 
conducted  for  this  thesis.  The  results  were  used  to  analyze 
the  two  objectives  which  served  as  guidance  for  this 
research.  In  addition,  the  results  obtained  were  used  to 
support  or  refute  the  research  hypothesis.  The  research 
objectives  will  be  restated,  and  results  and  analysis 
presented  in  the  remainder  of  this  chapter. 


Research  Objectives 

Research  Objective  1.  Validate  the  current 
forecasting  method  used  for  computing  tonnage 
estimates  to  derive  SDT  budget  requests. 

This  research  objective  was  accomplished  by  acquiring 
the  necessary  flying  hour  and  SDT  tonnage  data  from 
AFLC/DSXR  for  'he  40  quarters  from  FY  78/1  through  FY  87/4. 
As  explained  earlier,  these  data  were  for  SDT  tonnage 
shipments  by  MAC  and  MSC  to  USAFE  and  PACAF ;  therefore,  four 
distinct  data  sets  were  compiled.  In  addition,  the  actual 
AFLc/DSXR  SDT  tonnage  forecasts  for  the  six  quarters  from 
FY  88/1  through  FY  89/2  for  these  same  sets  were  obtained  to 
allow  for  accuracy  measurement  comparisons  between  the 
method  currently  used  and  the  models  developed  in  this 


research.  Finally,  data  concerning  the  number  of  quarters 


used  by  AFLC/DSXR  in  each  regression  forecast  for  each  data 


set  were  gathered  to  determine  if  the  number  of  quarters 
were  in  the  unstable  range  of  the  flying  hour  confidence 
intervals . 

Another,  and  final,  step  in  accomplishing  this  research 
objective  was  to  determine  if  the  flying  hour  parameter, 
in  AFLC/DSXR's  current  model  changed  statistically  after 
each  regression  iteration  was  conducted.  If  this  change  did 
occur,  the  model  was  considered  statistically  invalid.  As 
stated  earlier,  a  3-step  process  was  conducted  on  each  data 
set  to  determine  if  any  change  of  the  flying  hour  parameter 
occurred . 

1.  The  for  each  iteration  was  computed  using 
QUATTRO . 

2.  The  standard  error  (sp^)  for  each  was  computed 
using  QUATTRO. 

3.  Using  these  standard  errors,  95  percent  confidence 
intervals  for  for  each  iteration  were  established. 

Data  Analysis .  The  data  for  all  four  sets  were 
compiled  and  graphed  by  two  methods.  First,  SDT  tonnage, 
the  dependent  variable,  was  graphed  against  the  independent 
variable,  regional  flying  hours.  Secondly,  SDT  tonnage  was 
graphed  over  time. 

As  seen  in  Figure  2,  MAC  SDT  tonnage  shipped  to  USAFE, 
when  graphed  versus  flying  hours,  is  widely  dispersed  and 
appears  curvilinear. 
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Figure  2.  Regional  Flying  Hours  versus 
MAC  SDT  Tonnage  for  USAFE 


Figure  3  depicts  MAC  SDT  tonnage  shipped  to  USAFE  as  a 
time  series.  This  relationship  does  not  appear  to  be  linear 
and  constant,  but  rather  curvilinear. 

l'he  relationship  of  MSC  SDT  tonnage  shipped  to  USAFE  to 
flying  hours  is  depicted  in  Figure  4.  This  graph  suggests  a 
curvilinear  heteroscedas tic  relationship  of  tonnage  to 
flying  hours. 

Figure  5  exhibits  MSC  SDT  tonnage  shipped  to  USAFE  as  a 
time  series.  As  evidenced  by  the  graph,  the  series  is 
curvilinear  and  nonconstant.  There  is  no  indication  of  a 
constant  or  linear  relationship. 

Figure  6  graphically  relates  MAC  SDT  tonnage  shipped  to 
PACAF  to  regional  flying  hours.  As  the  figure  suggests, 
no  semblance  of  linearity  exists  between  the  two  variables. 
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Figure  3.  Quarterly  SDT  Tonnage  Shipped 
to  USAFE  by  MAC 
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Figure  4.  Regional  Flying  Hours  versus 
MSC  SDT  Tonnage  for  USAFE 
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Figure  5.  Quarterly  SDT  Tonnage  Shipped 
to  USAFE  by  MSC 


MAC  SDT  tonnage  shipped  to  PACAF  is  expressed  as  a  time 
series  in  Figure  7.  The  series  appears  linear  with  a  slight 
upward  trend. 


Figure  8  shows  MSC  SDT  tonnage  shipped  to  PACAF  plotted 
against  flying  hours.  Although  the  graph  appears  linear  at 
the  lower  range  of  flying  hours,  the  variance  within  the 
data  increases  with  the  amount  of  flying  hours. 

MSC  SDT  tonnage  shipped  to  PACAF  is  expressed  as  a  time 
series  in  Figure  9.  This  series,  although  linear  with  a 
slight  upward  trend  initially,  becomes  erratic  with  large 
variances  in  later  periods. 
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Figure  8.  Regional  Flying  Hours  versus 
MSC  SDT  Tonnage  for  PACAF 
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Figure  9. 


Quarterly  SDT  Tonnage  Shipped 
to  PACAF  by  MSC 


The  number  of  quarters  used  by  AFLC/DSXR  in  each  data 
set  regression  varied.  The  regressions  for  tiSC  SDT  tonnage 
shipped  to  both  USAFE  and  PACAF  used  40  quarters.  The 
regression  for  MAC  SDT  tonnage  shipped  to  PACAF  used  23 
quarters,  and  20  quarters  were  used  for  the  MAC  SDT  tonnage 
shipped  to  USAFE  regression. 

Equations  (26)  through  (29)  present  the  linear 
regression  models  used  by  AFLC/DSXP.  to  forecast  the  SDT 
tonnage  for  these  data  sets. 

The  linear  regression  used  for  MAC  SDT  tonnage  shipped 
to  USAFE  is  shown  in  equation  (26). 

Y  =  2332.137  +  0.093162X  !26) 

where;  Y  =  the  forecasted  SDT  tonnage 

X  =  the  flyinq  hours  associated  with  the  quarter 
being  forecasted 
n  =  20 
R  =  0.7175. 

Equation  (27)  expresses  the  linear  model  used  to 
forecast  the  MSC  SDT  shipped  to  USAFE. 

Y  =  -15,202.7  +  1  .  127496  X  •.2") 

where :  n  =  40 

R  =  0.5230. 

The  linear  regression  model  used  to  forecast  MAC  SDT 
tonnage  shipped  to  PACAF  is  shown  in  equation  (23) . 

Y  =  4374.293  +  0.0297  4  8X  t.23! 

where :  n  =  2  3 

R  =  0.1666. 
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Equation  (29)  was  used  as  the  linear  model  to  predict  MSC 
SDT  tonnage  shipped  to  PACAF . 


Y  =  -11 , 070  +  1 . 302432X 


(29) 


where:  n  =  40 

R  =  0.5235. 


The  forecasted  values  for  each  of  these  linear 
regressions  are  shown  in  Table  4. 


Table  4 

AFLC/DSXR  Forecasted'  Values  for  SDT  Tonnage 
MAC  SDT  MSC _ SDT 


FY 

USAFE 

PACAF 

USAFE 

PACAF 

88/1 

9,232 

6 , 141 

67 , 703 

44  ,  393 

88/2 

9 , 086 

6 , 123 

65 , 931 

43,583 

88/3 

9,773 

6 ,130 

74 , 243 

43 , 901 

88/4 

10 , 149 

6 , 131 

78,792 

43 , 964 

89/1 

9,609 

6,149 

72  ,  258 

44 . 747 

89/2 

9,611 

6 , 168 

72 , 288 

45  ,  562 

Curr e n t  Mode 1 

V a lid  a  t i on 

As  stated 

e  a  r  1  i 

er  in 

this  chapter,  a  3-step 

process 

was 

conducted  on 

each 

data 

set  to 

statistically  determine 

i  f 

changes  occurred  in  the 

flying  hour  parameter  after  each  regression  iteration.  The 
flying  hour  parameter  and  upper  and  lower  confidence 


interval  values  for  each  regression  iteration  within  each 
data  set  were  graphed.  These  values  are  tabulated  in 
Appendix  B.  The  graphs  indicate  if  any  of  the  confidence 
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intervals  within  a  particular  data  set  regression  failed  to 
overlap,  thus  invalidating  the  model. 

As  stated  earlier  in  the  thesis,  the  following 
hypothesis  test  was  used  to  determine  statistically  whether 
or  not  the  flying  hour  parameters  for  each  regression 
iteration  within  a  data  set  were  equal. 


H0  •* 

Pi, 

8  =  Pi , 9  = 

■  -  =  Pi , 40 

ha; 

Pi, 

,  i  *  Pi , 8  =  • 

•  •  =  P 1 ,  n 

where:  i  =  any  one  iteration  conducted  with  8 
to  40  periods  of  data 
n  =  total  number  of  iterations 
conducted  excluding  i 

Test  Statistic:  ±  t<z/2?(31 

where:  ta/2  =  the  value  of  the  test-statistic  for 

a  =  .05  and  n  -  2  degrees  of 
freedom  ( df ) . 

Rejection  Region:  Reject  Hq  if  any  two  of  the 

regression  iteration  confidence 
intervals  did  not  overlap. 

A  failure  of  any  two  confidence  intervals  to  overlap 
resulted  in  the  failure  to  support  t"he  null  hypothesis  (Hq) 
of  model  stability.  A  failure  of  this  type  statistically 
invalidated  AFLC/DSXR’s  current  model  as  a  forecasting  tool. 

Graphically,  a  failure  of  overlap  occurred  whenever  the 
lower  limit  of  one  confidence  interval  within  a  data  set  was 
la _ger  (horizontally  higher  on  the  graph)  than  the  upper 
limit  of  another  confidence  interval  within  the  same  data. 

Figure  10  shows  the  flying  hour  parameter  confidence 
intervals  for  Ml C  SDT  tonnage  shipped  to  USAFE. 
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Figure  10.  Confidence  Intervals  for  the 
Flying  Hour  Parameter  (MAC  to  USAFE) 


Based  on  the  graph  in  Figure  10,  H q  cannot  be  rejected  for 
MAC  SDT  tonnage  forecast  for  USAFE.  The  regression 
iteration  using  27  quarters  (lower  limit  of  0.0508)  did, 
however,  almost  fail  to  overlap  with  the  40-quarter 
regression  iteration  (upper  limit  of  0.0598'. 

Figure  11  depicts  the  parameter  confidence  intervals 
for  MSC  SDT  tonnage  forecasted  for  USAFE.  The  graph  in 
Figure  11  clearly  exhibits  several  confidence  intervals  that 
do  not  overlap.  Using  this  information,  Hq  was  rejected 
thus  resulting  in  the  conclusion  that  the  flying  hour 
parameter  changed  during  the  regression  iterations. 

Figure  12  shows  all  parameter  confidence  intervals  for 
MAC  SDT  tonnage  forecasted  for  PACAF  overlapping. 
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Quarters  in  the  Regression 


Figure  11.  Confidence  Intervals  for  the 
Flying  Hour  Parameter  (MSC  to  USAFE) 


Figure  12.  Confidence  Intervals  for  the 
Flying  Hour  Parameter  (MAC  to  PACAF) 


Figure  13  exhibits  parameter  confidence  mcervdls  for 
MSC  SDT  tonnage  forecast  for  PACAF.  The  graph  shows,  again, 
several  confidence  intervals  not  overlapping.  These 
failures  result  in  a  rejection  of  Ho  leading  to  the 
conclusion  that  the  parameter  changed  during  the  regression 
iterations . 


40  30  20  1 0 

35  .25  .15 

Quarters  in  the  Regression 


Lo  Limit 
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Uc  Limit 


Figure  13.  Confidence  Intervals  for  the 
Flying  Hour  Parameter  (MSC  to  PACAF) 


In  summary,  the  statistical  test  employed  to  determine 
the  validity  of  the  current  model  resulted  in  the  two  MSC 
SDT  data  sets,  in  fact,  being  statistically  invalidated; 
however,  the  same  test  failed  to  statistically  invalidate 
the  two  MAC  SDT  data  sets.  This  result  does,  however,  imply 
the  current  model,  in  general,  is  perhaps  an  ineffective 
forecasting  tool. 
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Research  Objective  2.  If  the  curreni  method's 
validity  was  not  supported,  develop  a  new 
forecasting  model,  using  the  same  input  data,  that 
would  produce  more  accurate  and  reliable  tonnage 
estimates . 

As  noted  above,  the  validity  of  the  current  method  is 
suspect.  Therefore,  the  second  research  objective  was 
accomplished  in  two  steps.  First,  Box-Jenkins  (BJ)  time 
series  forecasting  models  were  built  for  each  of  the  four 
data  sets.  Secondly,  the  accuracy  of  these  models  was 
evaluated  using  the  mean  absolute  percent  error  (MAPE) .  The 
MAPE  was  used  because  it  enabled  the  size  of  the  error  in  a 
particular  period  to  be  related  to  the  actual  tonnage  in  the 
same  period;  thus  the  MAPE  offers  the  user  the  error  as  a 
size  percentage  of  the  actual  tonnage. 

As  stated  earlier  in  the  methodology  chapter,  BJ  time 
series  forecasting  was  used  to  develop  new  models  for  each 
data  set  to  forecast  SDT  tonnage.  This  choice  was  made 
because  BJ  enables  patterns  to  be  identified  in  the  history 
of  a  time  series  and  uses  the  patterns  to  build  the 
appropriate  model. 

As  stated  earlier  in  Chapter  II,  four  distinct  steps 
were  followed  in  building  these  BJ  models. 

1.  Identification  of  any  patterns  in  the  time  series. 

2.  Model  specification  based  on  these  identified 
patterns . 
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3.  Diagnostic  tests  to  ensure  the  appropriate  model  is 
specified . 

4.  Hypothesis  testing  and  forecasting. 

Pattern  Identification.  Data  plots,  auto¬ 
correlation  functions  (ACFs),  and  partial  autocorrelation 
functions  (PACFs)  were  graphed  and  analyzed  for  each  of  the 
time  series  data  sets  to  identify  any  underlying  patterns 
within  the  series.  The  araphs  shown  in  Figures  14  through 
17  repeat  the  time  series  graphs  for  MAC  and  MSC  SDT  tonnage 
shipped  to  both  USAFE  and  PACAF.  Using  these  graphs  and  the 
ACFs  and  PACFs  computed  with  the  TIMES  forecasting  software 
package,  any  underlying  patterns  in  the  history  of  the  time 
series  were  identified.  The  ACFs  and  PACFs  are  located  in 
Appendix  C. 
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Figure  14.  Quarterly  SDT  Tonnage  Shipped 
to  USAFE  by  MAC 
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Figure  15.  Quarterly  SDT  Tonnage  Shipped 
to  USAFE  by  M SC 
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76/1  80/1  82/1  84/1  86/1 
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Fiscal  CX/arter 

Figure  16.  Quarterly  SDT  Tonnage  Shipped 
to  PACAF  by  MAC 
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73/1  80/1  82/1  84/1  86/1 

79/1  81  /I  83/1  85/1  87/1 

Fiscal  Quarter 

Figure  17.  Quarterly  SDT  Tonnage  Shipped 
to  PACAF  by  MSC 

All  the  graphs,  with  the  exception  of  the  MAC  SDT 
tonnage  shipped  to  PACAF,  exhibit  some  form  of 
nonstationarity .  This  nonstationarity  implied  that  when 
constructing  the  models  some  degree  of  differencing  was 
necessary  to  make  the  data  stationary. 

In  referring  to  Appendix  C,  the  ACF  and  PACF  for  MAC 
SDT  tonnage  shipped  to  USAFE  indicated  a  need  to  difference 
the  data.  The  degree  of  differencing  required  was 
determined  by  comparing  the  value  "to  test  whether  this 
series  is  white  noise"  (sometimes  referred  to  as  the 
Q-statistic)  for  the  original  series  to  the  two  differenced 
series.  Since  the  value  for  the  original  series,  44.98,  was 
larger  than  the  values  for  the  other  two  series  (44.14  and 
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35.63),  differencing  of  order  two  was  conducted.  The  ACF 
and  PACF  for  the  two  difference  series  were  used  to 
determine  what  parameters  should  be  included  in  the  model . 
Since  significant  spikes  were  apparent  in  both  functions,  AR 
and  MA  parameters  were  included. 

The  ACF  and  PACF  for  MSC  SDT  tonnage  shipped  to  USAFE 
indicated  a  need  to  difference.  In  comparing  the  Q- 
statistic  for  the  original  series  to  the  two  differenced 
series,  the  value  for  the  original  series,  177.28,  was 
larger  than  both  the  values  for  the  remaining  two  series 
(28.78  for  one  difference  and  58.16  for  two  differences). 
Since  the  r>ne  difference  series  possessed  the  smallest  Q- 
statistic,  differencing  of  order  one  was  conducted  in 
building  the  model  for  this  data  set.  The  ACF  and  PACF  for 
the  one  difference  series  were  used  to  determine  what 
parameters  should  be  included  in  the  model.  Since 
significant  spikes  were  evident  in  both  functions  at  the 
first  lag,  AR  and  MA  parameters  were  included. 

The  ACF  and  PACF  for  MAC  SDT  tonnage  shipped  to  PACAF 
indicated  no  need  to  difference.  The  Q-statistic  for  the 
original  series,  15.86,  was  smaller  than  the  values  for  the 
remaining  two  series  (15.99  and  19.27) .  The  ACF  and  PACF 
for  the  original  series  were  used  to  determine  what 
parameters  should  be  included  in  the  model.  Since  there  was 
a  significant  spike  in  the  ACF  and  the  PACF  tailed  off,  only 
an  MA  parameter  was  included. 
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The  ACF  and  PACF  for  MSC  SDT  tonnage  shipped  to  PACAF 
indicated  a  need  to  difference.  The  Q-statistic  for  the 
original,  113.48,  was  larger  than  the  values  for  the 
remaining  two  series  (19.98  and  44.49).  The  ACF  and  PACF 
for  the  one  difference  series  were  used  to  determine  what 
parameters  should  be  included  in  the  model.  Since  there 
were  significant  spikes  in  the  ACF  and  the  PACF,  AR  and  MA 
parameters  were  included. 

Model  Specif ication .  Using  the  information 
gathered  in  the  pattern  identification  phase,  appropriate 
models  were  specified.  Again,  keeping  in  mind  that  the  BJ 
model  building  process  is  iterative,  Equations  (30)  through 
(33)  present  the  final  models  selected. 

Equation  (30)  shows  the  ARIMA  (1,2,2)  model  used  for 
MAC  SDT  tonnage  shipped  to  USAFE. 

Yt  =  Yt-2  +  *2<Yt-2  ~  Yt- 4 )  “  ©let-l 

-  ©2et-2  +  et  (3°- 

where:  Yt  =  the  forecasted  time  series  values 

Yt-i  =  historical  value  of  the  time  series 

=  the  AR  parameter  at  period  i 

0j_  =  the  MA  parameter  at  period  i 

et-i  =  the  error  associated  with  period  t  -  i. 

This  model  accounted  for  the  need  to  difference  the  original 
data  set  twice.  The  model  also  used  the  patterns  associated 
with  the  ACF  and  PACF.  The  MA  parameters  included  in  the 
model  were  deemed  necessary  due  to  the  two  signifi  nt 
spikes  at  lags  one  and  two  in  the  ACF,  and  the  AR  parameter 
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was  included  to  account  for  the  significant  spike  at  lag 
three  m  the  PACF . 

Equation  (31)  presents  the  ARIMA  (1,1,1.  model 
associated  with  the  MSC  SDT  tonnage  shipped  to  USAFE. 


Yt  =  Yt-i  +  ♦itYt-!  -  Yt-2)  -  ©iet-i  +  et  (31) 

where:  =  the  forecasted  time  series  values 

Y^-i  =  historical  value  of  the  time  series 

4>i  =  the  AR  parameter  at  period  i 

©i  =  the  MA  parameter  at  period  i 

e^-i  =  the  error  associated  with  period  t  -  i. 


This  model  accounted  for  the  need  to  difference  the  original 
data  once.  It  also  utilized  the  appropriate  parameters  to 
coincide  with  the  ACF  and  PACF.  The  MA  parameter  was 
included  to  account  for  the  maior  significant  spike  at  laa 
one  in  the  ACF.  The  AR  parameter  was  used  to  explain  the 
large  spike  at  lag  one  in  the  PACF. 

The  ARIMA  (0,0,1)  model  in  equation  (32)  was  used  with 
the  data  for  MAC  SDT  tonnage  shipped  to  PACAF. 

Y^  =  e^  -  +  U  (32) 

where:  Y t  =  the  forecasted  time  series  values 
©i  =  the  MA  parameter  at  period  l 
e(--i  =  the  error  associated  with  period  t  -  i 
U  =  the  mean  of  the  time  series. 

The  model  in  equation  (32)  accounted  for  the  necessary 
patterns  identified  in  the  data  plot,  ACF,  and  PACF.  Only 
the  one  MA  parameter  was  included  to  account  for  the 
significant  spike  at  lag  one  in  the  ACF  and  the  tailing -off 


of  the  PACF. 


The  ARIMA  (1,1,2)  model  in  equation  (33)  was  used  for 
M SC  SDT  tonnage  shipped  to  PACAF. 

Yt  =  Yt-l  +  *l(Yt-l  -  Yt-2) 

"  03et-3  “  ©4et-4  +  et  (33) 

where:  Y^  =  the  forecasted  time  series  values 

Yt-i  =  historical  value  of  the  time  series 

=  the  AR  parameter  at  period  i 

©i  =  the  MA  parameter  at  period  i 

et-i  =  the  error  associated  with  period  t  -  i. 

This  model  accounted  for  the  need  to  difference  and  the 
necessary  patterns  identified  previously  in  the  time  series. 
The  two  MA  parameters  were  used  to  explain  the  spikes  at 
lags  three  and  four  in  the  ACF,  while  the  AR  parameter  was 
used  to  account  for  the  spike  at  lag  one  in  the  PACF. 

The  estimated  parameters  and  models  are  in  Equations 
(34)  through  (37)  located  in  the  New  Model  Forecasts  section 
below . 

Diagnostic  Testing.  Eight  diagnostic  tests  were 
conducted  on  each  of  the  four  specified  models.  These  tests 
wf>re  conducted  using  TIMES.  The  necessary  output  from  TIMES 
concerning  these  tests  is  located  in  Appendix  D.  As  stated 
earlier  in  the  methodology  chapter,  the  BJ  model  building 
process  is  iterative.  Therefore,  more  than  one  model  for 
each  data  set  was  specified.  The  results  of  the  diagnostic 
tests  for  each  of  the  models  for  each  data  set  were  compared 
to  each  other.  Only  the  diagnostic  results  from  the  models 
chosen  for  each  data  set  are  discussed  below. 
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The  diagnostic  output  for  MAC  SDT  tonnage  to  USAFE 
showed  the  following: 

1.  The  mean  of  the  residuals  was  47.36  and  the 
standard  error  remained  constant  compared  to  the  other 
models  that  were  specified  for  this  time  series. 

2.  The  residual  ACF  showed  only  a  small  "spike"  at  the 
third  lag,  indicating  that  not  all  of  the  disturbance  in  the 
time  series  was  accounted  for  by  the  specified  model. 

3.  The  Q-statistic  was  3.276  and  when  compared  with  a 
chi-square  variable  with  21  degrees  of  freedom,  the 
associated  alpha  (a)  value  fell  between  99  and  99.5  percent. 
Thus,  it  was  concluded  with  99  to  99.5  percent  confidence, 
th_.c  the  residual  autocorrelations  were  not  significant 
implying  that  the  model  accounted  for  all  of  the  significant 
disturbance  in  the  time  series 

4.  The  cumulative  periodogram  of  the  residuals  was 
approximately  linear  and  remained  within  the  expected 
probability  limits,  thus  implying  that  the  residuals  were 
normalized  about  a  mean  of  zero. 

5.  The  histogram  of  the  residuals  showed  a  mound- 
shaped  graph  centered  about  zero,  implying  che  residuals  met 
the  standard  assumption  of  the  error  term  --  being 
normalized  and  centered  about  a  mean  of  zero. 

6.  The  3 S E  was  13,11  3,000  and  the  MSE  was  ^''",980. 
These  values  when  compared  to  similar  33E  and  MSE  values  t 
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other  models  specified  for  this  data  set  at  33  degrees  of 
freedom  were  smaller  and  therefore  more  acceptable. 

7.  The  Power  Spectrum  was  horizontal  and  comparatively 
smooth.  This  finding  indicated  no  problems  with  the  error 
term  in  the  model. 

8.  The  Schwarz  statistic,  or  Bayesian  Information 
Criterion  (BIC),  was  equal  to  17,318,826  and  was  smaller 
than  the  BIC  for  the  other  models  specified  for  this  time 
series,  implying  that  the  model  selected  was  the  best 

f i tting . 

The  diagnostic  output  for  HSC  SDT  tonnage  to  USAFE 
showed  the  following: 

1.  The  mean  of  the  residuals  was  834.31  and  the 
standard  error  remained  constant  compared  to  the  other 
models  that  were  specified  for  this  time  series.  The  mean 
was  large,  but  acceptable  based  on  the  original  time  series 
being  in  the  tens  of  thousands. 

2.  The  residual  ACF  showed  no  spikes,  indicating  the 
model  specified  accounted  for  all  the  patterns  in  the  time 
series  . 

3.  The  Q-statistic  was  13.505  and  when  compared  wi.h  a 
chi-square  variable  with  22  degrees  of  freedom,  the 
associated  alpha  (a)  value  fell  between  90  and  95  percent 
Thus,  it  was  concluded  with  90  to  95  percent  confidence, 
that  the  residual  autocorrela*'  ions  were  not  significant 


implying  that  the  model  accounted  for  all  of  the  significant 
disturbance  in  the  time  series. 

4.  The  cumulative  periodogram  of  the  residuals  was 
approximately  linear  and  remained  within  the  expected 
probability  limits,  thus  implying  that  the  residuals  were 
normalized  about  a  mean  of  zero. 

5.  The  histogram  of  the  residuals  showed  a  mound¬ 
shaped  graph  centered  about  zero,  implying  the  residuals  met 
the  standard  assumption  of  the  error  term  --  being 
normalized  and  centered  about  a  mean  of  zero. 

6.  The  S2E  was  2,591,500,000  and  the  MSE  was 
71,937,000.  These  values  when  compared  to  similar  SSE  and 
MSE  values  of  other  models  specified  for  this  data  set  at  36 
degrees  of  freedom  were  smaller  and  therefore  more 
acceptable . 

7.  The  Power  Spectrum  was  horizontal  and  relatively 
smooth.  This  finding  indicated  no  problems  with  the  error 
term  in  the  model . 

8.  The  BIC  was  equal  to  3,116,407,298  and  was  smaller 
than  the  BIC  for  the  other  models  specified  for  this  time- 
series,  implying  that  the  model  selected  was  the  best 
fitting. 

The  diagnostic  output  for  MAC  SDT  tonnage  to  PACAF 
showed  the  following; 

1.  The  mean  of:'  the  residuals  was  6.5003  and  the 
standard  error  remained  constant  compared  to  the  other 


models  that  were  specified  for  this  time  series. 

2.  The  residual  ACF  showed  no  jpikes,  indicating  the 
model  specified  accounted  for  all  the  patterns  in  the  time 
series . 

3.  The  Q-statistic  was  9.4246  and  when  compared  with  a 
chi-square  variable  with  22  degrees  of  freedom,  the 
associated  alpha  (a)  value  fell  between  99  and  99.5  percent. 
Thus,  it  was  concluded  with  99  to  99.5  percent  confidence, 
that  the  residual  autocorrelations  were  not  significant 
implying  that  the  model  accounted  for  all  of  the  significant 
disturbance  in  the  time  series. 

4.  The  cumulative  periodogram  of  the  residuals  was 
linear  and  remained  within  the  expected  probability  limits, 
thus  implying  that  the  residuals  were  normalized  about  a 
mean  of  zero. 

5.  The  histogram  of  the  residuals  showed  a  mound¬ 
shaped  graph  centered  about  zero,  implying  the  residuals  met 
the  standard  assumption  of  the  error  term  --  being 
normalized  and  centered  about  a  mean  of  zero. 

6.  The  SSE  was  6,011,300  the  MSE  was  153,190. 

These  values  when  compared  to  similar  SSE  and  MSE  values  cf 
other  models  specified  for  this  data  set  at  33  degrees  of 
freedom  were  smaller  and  therefore  more  acceptable. 

7.  The  Power  Spectrum  was  horizontal  and  smooth.  This 
finding  indicated  no  problems  with  the  error  term  in  the 
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model  . 


8.  The  BIC  waj  equal  to  6,592,041  and  was  smaller  than 
the  BIC  for  the  other  models  specified  for  this  time  series, 
implying  that  the  model  selected  was  the  best  fitting. 

The  diagnostic  output  for  MSC  SDT  tonnage  to  PACAF 
showed  the  following: 

1.  The  mean  of  the  residuals  was  172.20  and  the 
standard  error  remained  constant  compared  to  the  other 
models  that  were  specified  for  this  time  series.  The  mean 
was  large,  but  acceptable  based  on  the  original  series  being 
in  the  tens  of  thousands. 

2.  The  residual  ACF  showed  no  spikes,  indicating  the 
model  specified  accounted  for  all  the  patterns  in  the  time 
series  . 

3.  The  Q-statistic  was  11.469  and  when  compared  with  a 
chi-square  variable  with  21  degrees  of  freedom,  the 
associated  alpha  (a)  value  fell  between  95  and  97.5  percent. 
Thus,  it  was  concluded  with  95  to  97.5  percent  confidence, 
that  the  resilual  autocorrelations  were  not  significant 
implying  that  the  model  accounted  for  all  of  the  significant 
disturbance  in  the  time  series. 

4.  The  cumulative  periodogram  of  the  residuals  was 
approximately  linear  and  remained  within  the  expected 
probability  limits,  thus  implying  that  the  residuals  were 
normalized  about  a  mean  of  zero. 

5.  The  histogram  cf  the  residuals  showed  a  mouni- 
shaped  graph  centered  about  zero,  implying  the  residuals  mef 
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the  standard  assumption  of  the  error  term  --  being 
normalized  and  centered  about  a  mean  of  zero. 

6.  The  SSE  was  735,760,000  and  the  MSE  was  21,022,000. 
These  values  when  compared  to  similar  SSE  and  MSE  values  of 
other  models  specified  for  this  data  set  at  35  degrees  of 
freedom  were  smaller  and  therefore  more  acceptable. 

7.  The  Power  Spectrum  was  horizontal  and  smooth.  This 
finding  indicated  no  problems  with  the  error  term  in  the 
model . 

8.  The  BIC  was  equal  to  970,265,687  and  was  smaller 
than  the  BIC  for  the  other  models  specified  for  this  time 
series,  implying  that  the  model  selected  was  the  best 
fitting. 

These  diagnostic  tests  confirmed  that  the  BJ  models 
chosen  for  forecasting  the  four  time  series  were 
statistically  valid  and  the  most  acceptable  models  of  those 
specified . 

New  Model  Forecasts.  After  the  specified  models 
passed  the  diagnostic  tests  with  stability,  forecasts  were 
conducLc'3  for  FY  88/1  through  FY  89/2. 

The  ARIMA  (1,2,2)  model  used  to  forecast  MAC  SDT 
tonnage  to  USAFE  is  expressed  in  equation  (34).  The 
equation  includes  the  numerical  parameters. 

Y t  =  Yt-2  -  0 . 50572 (Yt_2  -  Yt_4) 

-  1.4292et_1  +  0.73174et_2  +  et  (34) 
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The  ARIMA  (1,1,1)  model  including  the  numerical 
parameters  used  to  forecast  MSC  SDT  tonnage  to  USAFE  is 
shown  in  equation  (35). 

Yt  =  Yt_i  -  0.809(Yt_i  -  Yt_2)  +  0.5753et_i  +  et  (35) 

MAC  SDT  tonnage  shipped  to  PACAF  was  forecasted  using 
the  ARIMA  (0,0,1)  model  expressed  in  equation  (36). 

Y{-  =  e^  +  0.40823et-l  +  6064.4  (3b) 

The  ARIMA  (1,1,2)  model  shown  in  equation  (37)  was  used 
to  forecast  MSC  SDT  tonnage  to  PACAF. 

Yt  =  Yt-1  ~  0.26403 (Yt_!  -  Yt_2) 

-  0 . 074567et_3  +  0.1363bet_4  =  et  (37) 

The  forecasted  values  for  each  of  these  time  series  are 
shown  in  Table  5.  These  can  be  compared  with  the  AFLC/DSXR 
forecasts  in  Table  4  and  the  actual  tonnage  in  Appendix  A. 

Table  5 

Box-Jenkins  Forecasted  Values  for  SDT  Tonnage 


MAC  S  DT  MSC SDT 


FY 

USAFE 

PACAF 

USAFc, 

PACAF 

88/1 

8  ,  516 

5,909 

70,021 

37 , 004 

88/2 

8 ,631 

6  064 

69 , 589 

37,380 

C3/3 

8  ,  890 

£>  ,  064 

69,939 

36 , 093 

88/4 

8  ,  542 

6 , 064 

69 ,656 

36,939 

89/1 

8 , 121 

6 , 064 

69,885 

36,717 

89/2 

8 , 007 

6 , 064 

69 , 700 

36,776 
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Accuracy  Measurement.  As  discussed  in  Chapter 
III,  the  mean  absolute  percent  error  (MAPE)  was  used  to 
measure  the  accuracy  of  AFLC/DSXR's  forecasting  model  and 
the  BJ  models  built  during  this  research.  The  forecast 
errors  for  AFLC/DSXR's  models  were  computed  by  subtracting 
the  actual  SDT  tonnage  from  the  forecasts.  Likewise,  the 
errors  for  the  BJ  models  developed  in  this  research  were 
computed  by  subtracting  the  actual  tonnage  from  the 
predictions.  Table  6  includes  the  MAPE  for  each  of  the 
models . 


Table  6 


MAPE 

Evaluations  of  the 

Forecast  Models 

MAC 

SDT 

MSC_ 

SDT 

Model 

USAFE 

PACAF 

USAFE 

PACAF 

AFLC/DSXR 

32.95 

35.40 

6 .61 

7 .87 

AFLC/DSXR 

3 , 037 

2,151 

4 , 295 

3 , 0  8  L 

(in  tons) 

Box- Jenkins 

19.31 

33 . 25 

6 . 97 

18.85 

Box- Jenkins 

1 ,780 

2 , 021 

4  ,  529 

7  ,  339 

(in  tons) 

The  tonnage  values  in  Table  6  are  the  absolute  error  in 
tonnage.  In  effect,  they  are  the  percentage  of  the  average 
tonnage  shipped  for  that  particular  data  set. 


73 


Summary 

This  chapter  presented  the  results  and  analysis  of  this 
research.  These  results  allowed  for  the  findings  and 
implications,  included  in  Chapter  V,  to  be  drawn.  These 
findings  and  implications  are  based  solely  on  the  results 
and  analysis  of  this  research. 
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V. 


Findings  and  Impli cations 

This  chapter  is  devoted  to  the  discussion  of  the 
findings  resulting  from  this  research  and  the  implications 
of  these  findings  on  the  forecasting  of  AFL.C  SDT  tonnage. 

F indings 

Research  Objective  1.  Validate  the  current 
forecasting  method  used  for  computing  tonnage 
estimates  to  derive  SDT  budget  requests. 

As  stated  in  Chapter  IV,  this  research  objective  was 
accomplished  by  data  analysis  and  a  3-step  model  validation 
process . 

Data  Analysis.  The  findings  from  graphically 
analyzing  the  four  data  sets  were  that  none  of  these  sets 
were  entirely  appropriate  for  use  with  linear  regression. 
Graphing  the  SDT  tonnage  versus  flying  hours  showed  all  four 
sets  were  curvilinear,  cone-shaped,  or  blocked,  and 
therefore  inappropriate  for  use  with  linear  regression.  In 
addition,  when  SDT  tonnage  was  graphed  as  a  time  series,  the 
same  findings  were  reached.  The  data  in  these  four  sets  did 
not  lend  themselves  for  use  with  linear  regression. 

Current  Model  Validation.  The  findings  from  the 
3-step  model  validation  process  were  varied.  All  flying 
hour  parameter  confidence  intervals  for  MAC  SDT  tonnage 
shipped  to  both  USAFE  and  PACAF  overlapped.  There  were, 
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however,  several  parameter  confidence  intervals  failed  to 
overlap  for  MSC  SDT  tonnage  shipped  to  USAFE  and  PACAF. 

Based  on  these  failures,  the  null  hypothesis  that  the 
parameters  were  stable  after  each  regression  was  rejected. 
The  iterative  regression  process  utilized  by  AFLC/DSXR  to 
predict  MSC  SDT  tonnage  shipped  to  both  USAFE  and  PACAF  was 
therefore  statistically  invalidated  based  on  a  95  percent 
confidence  level.  This  means  that  Hq  was  rejected  as  being 
false  with  only  a  five  percent  chance  of  rejecting  Hq  and  it 
being  true.  This  research,  however,  resulted  in  a  failure 
to  reject  the  null  hypothesis  for  predictinq  MAC  SDT  tonnage 
shipped  to  USAFE  and  PACAF. 

The  varying  number  of  quarters  used  in  each  regression 
also  contributed  to  the  instability  of  the  model.  A! thouch 
the  regressions  for  MSC  SDT  tonnage  shipped  to  USAFE  and 
PACAF  used  40  quarters  and  were  therefore  in  the  more  stable 
range  of  the  confidence  interval  graphs,  the  regressions  for 
MAC  SDT  tonnage  shipped  to  USAFE  and  PACAF  fell  in  the  more 
unstable  ranges.  The  regressions  for  MAC  SDT  tonnage 
shipped  to  USAFE  and  PACAF  used  20  and  23  quarters 
respectively.  These  two  regressions  fell  in  the  conf  idem.'-3 
interval  ranges  unere  the  intervals  began  increasing  in  size 
to  account  for  the  same  probability  of  the  parameter  values 
falling  within  the  interval  limits.  In  addition,  it  was  in 
these  ranges  where  the  intervals  began  to  fluctuate  more 
widely  and  become  unstable.  This  increase  in  range  size 


SI 


implies  that  larger  degrees  of  error  are  necessary  in  order 
for  the  confidence  interval  to  include  the  actual  tonnage 
value,  based  on  a  95  percent  confidence  level. 

In  summary,  this  first  research  objective  was  met  in 
that  the  current  method  was  statistically  invalidated.  It 
is  evident,  based  on  u  95  percent  confidence  level,  that  the 
flying  hour  parameteis  for  the  MSC  SDT  data  sets  are  in  fact 
unstable  after  certain  regressions  used  during  AFLC/DSXR's 
iterative  forecasting  method.  This  finding  undoubtedly 
casts  suspicion  on  the  method,  in  general,  as  a  forecasting 
tool . 

Research  Objective  2.  If  the  current  method's 
validity  was  not  supported,  develop  a  new 
forecasting  model,  using  the  same  input  data,  that 
would  produce  more  accurate  and  reliable  tonnage 
estimates . 

This  research  bjective,  as  stated  earlier,  was 
accomplished  in  two  steps.  First,  Box-Jenkins  (BJ)  time 
series  forecasting  models  were  built  for  each  of  the  four 
data  sets.  Secondly,  the  accuracy  of  these  models  was 
evaluated  and  compared  to  that  of  the  current  model  using 
the  mean  absolute  percent  error  (MAPE) . 

Mode  1  Building.  In  building  BJ  models  for  each  of 
these  four  times  series,  the  '’-step  BJ  model  building 
process,  discussed  in  Chapter  III,  was  conducted.  The 
results  were  four  separate  ARIMA  models,  one  for  each  time 
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series.  These  ARIMA  models,  built  using  patterns  identified 
in  the  history  of  the  time  series,  were  then  used  to 
forecast  SDT  tonnage  for  FY  88/1  through  FY  89/2. 

Accuracy  Measurement .  The  forecasts  made  using 
these  newly  developed  BJ  models  were  compared  to  AFLC/DSXR's 
forecasts  using  the  MAPE.  The  results  of  this  comparison 
were  varied  and  are  included  previously  in  Table  6. 

The  MAPE  from  the  BJ  model  forecast  for  MAC  SDT  tonnage 
shipped  to  both  USAFE  and  PACAF  were  smaller  than  the  MAPE 
for  the  AFLC/DSXR  models.  These  were  the  two  data  sets 
where  the  AFLC/DSXR  method  could  not  be  statistically 
invalidated.  This  supports  the  fact  that  the  BJ  models 
developed  for  these  data  sets  provide  more  accurate 
forecasts  than  the  "valid"  AFLC/DSXR  models.  The  MAPE  from 
the  BJ  model  forecast  for  MSC  SDT  tonnage  shipped  to  USAFE 
and  PACAF  were  larger  than  those  for  the  AFLC/DSXR  models. 
The  BJ  models  for  these  two  data  sets  are,  however,  valid 
models,  whereas  the  AFLC/DSXR  method  used  for  the  same  data 
sets  was  statistically  invalidated. 

The  first  finding  above  met  the  research  objective  of 
developing  a  new  forecasting  model  that  was  a  more  accurate 
predictor  of  SDT  tonnage  than  the  current  model.  The  MAPE 
values,  however,  were  large  (19.31  percent  for  MAC  SDT 
tonnage  forecast  for  USAFE  and  33.25  percent  for  MAC  SDT 
tonnage  forecast  for  PACAF)  in  relation  to  the  mean  value  of 
each  corresponding  time  series.  These  values  imply  that  the 
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mean  absolute  error  of  the  forecast  is  19.31  percent  of  the 
mean  actual  tonnage  for  the  MAC  SDT  tonnage  shipped  to 
USAFE,  and  33.25  percent  of  the  mean  actual  tonnage  for  the 
MAC  SDT  tonnage  shipped  to  PACAF.  These  values  translate 
into  1,780  and  2,021  short  tons  for  the  respective  time 
series.  Yet,  even  with  these  relatively  large  MAPE  values, 
the  BJ  models  developed  are  more  accurate  in  predicting 
future  values  of  these  two  time  series  than  the  AFLC/DSXR 
models  which  yield  errors  of  3,037  and  2,151  short  tons 
respectively . 

One  final  note  here  is  that  the  AFLC/DSXR  models  for 
these  two  time  series  were  not  statistically  invalidated 
during  this  research  and  were  therefore  considered  "valid" 
models  based  on  the  95  percent  confidence  level  chosen  for 
this  research.  It  is  therefore  an  important  finding  that 
the  BJ  models  developed  for  these  time  series  were  more 
accurate  than  the  "valid"  AFLC/DSXR  models  for  the  same 
series . 

The  second  finding  noted  above  that  the  BJ  MAPE  values 
for  the  two  MSC  SDT  time  series  were  larger  than  the 
corresponding  MAPE  values  for  the  AFLC/DSXR  models  was 
attributed  to  the  "naivete"  of  the  BJ  time  series  modeling 
process.  Due  to  SDT  funding  cuts,  MAC  SDT  shipments  were 
diverted  to  MSC  SDT  shipments  causing  unexpected  increases 
in  the  MSC  SDT  shipments  between  FY  83/1  and  FY  89/2.  The 
BJ  process  uses  history  of  a  time  series  to  forecast  the 
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future  of  that  time  series  and  was  therefore  unable  to 
detect  or  predict  these  increases.  Even  though  the  BJ 
models  developed  for  these  two  MSC  SDT  time  series  were  less 
accurate  than  the  AFLC/DSXR  models  for  the  same  series  based 
on  the  MAPE  values,  they  were  both  relatively  accurate.  The 
BJ  MAPE  values  (6.97  for  MSC  SDT  tonnage  forecast  for  USAFE 
and  18.85  for  MSC  SDT  tonnage  forecast  for  PACAF)  imply  that 
the  mean  absolute  error  is  6.97  percent  of  the  mean  actual 
MSC  SDT  tonnage  shipped  to  USAFE,  and  18.35  percent  cf  the 
mean  actual  MSC  SDT  tonnage  shipped  to  PACAF.  These  values 
translate  into  4,529  and  7,389  measurement  tons  for  the 
respective  time  series.  The  values  are  relatively  small 
forecast  errors  considering  the  average  tonnage  shipped  tor 
the  40  quarters  in  the  time  series  was  64,980  and  39,198 
measurement  tons  respectively. 

One  final  note  here  is  that  although  the  AFLC/DSXR 
models  for  these  two  time  series  were  more  accurate  than  th* 
BJ  models  developed  for  these  same  series,  these  two 
AFLC/DSXR  models  were  the  two  that  were  statistically 
invalidated  during  this  research.  Therefore,  a  choice  must 
be  made  between  a  statistically  invalid  model  that  is 
marginally  more  accurate  or  a  statistically  valid  model  that 
is  marginally  less  accurate. 

In  summary,  this  second  research  objective  was  met  by 
developing  the  appropriate  BJ  models.  The  BJ  models  for  the 
two  i i AC  SDT  time  series  were  more  accurate  than  the 
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AFLC/DSXR  models  for  the  same  series  based  on  the  MAPE .  The 


BJ  models  developed  for  the  two  MSC  SDT  time  series, 
although  marginally  less  accurate  than  the  AFLC/DSXR  models 
for  the  same  series,  were  valid  models,  whereas  the 
AFLC/DSXR  models  were  statistically  invalidated  during  the 
course  of  this  research. 

Impl ications 

The  results  presented  here  offer  some  valuable  insight 
into  the  management  of  the  AFLC  SDT  budget  forecasting 
process . 

Continual  Analysis  and  Updating.  Although  the  SDT 
tonnage  data  sets  evaluated  in  this  research  may  have  been 
linear  in  the  early  1970's  when  AFLC/DSXR's  model  was  built, 
that  is  not  the  case  at  the  present  time.  Forecasting 
inputs  and  outputs  must  be  continually  analyzed  and  updated 
according  to  the  current  conditions  of  the  data  being 
forecasted.  Using  a  forecasting  method  inappropriate  for  a 
particular  data  set  serves  no  purpose  and  usually  provides 
inaccurate  results.  AFLC/DSXR  has  used  the  same  iterative 
linear  regression  model  for  the  past  several  years.  It  is 
apparent,  however,  by  reviewing  the  plots  of  the  data  in 
Chapter  IV  that  there  are  no  linear  relationships. 
AFLC/DSXR's  current  approach  is  therefore  inappropriate  for 
the  data  sets  it  forecasts. 

Accurate  Model  Need.  The  need  for  an  accurate  and 
statistically  valid  forecasting  model  to  predict  SDT  tonnage 
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is  necessary  as  a  result  of  AFLC/DSXR's  current  model  being 
statistically  invalidated  in  this  research.  Based  on  the 
fact  that  there  are  external  factors  affecting  these  time 
series  such  as  the  diversion  of  cargo  from  air  to  sea 
movement  and  the  delay  of  some  shipments,  self -projecting 
methods  (i.e.  Box-Jenkins)  become  less  effective. 
Econometric  models  are  more  apt  to  detect  these  external 
effects  and  therefore  become  better  predictors  of  such  time 
series  . 

Data  Base  Size  Increase.  Any  forecasting  method 
whether  linear  regression,  Box-Jenkins,  or  econometric 
requires  a  vast  amount,  usually  at  least  50  to  60  periods, 
of  historical  data.  The  amount  of  data  (40  quarters)  used 
in  this  research,  although  statistically  acceptable,  was 
limited.  Ideally,  50  to  60  periods  (for  this  research, 
quarters)  should  be  used  to  more  accurately  forecast  any 


given  time  series. 

This 

implication 

wou 

Id  r 

equire 

AFLC/DSXR 

to  retain  13  to  15 

years 

of  data  rather 

than 

ten . 

Alternatively,  monthly  da 

ta  for  four 

to 

five 

years 

could  be 

used.  If,  however,  an  econometric  model  was  employed, 
additional  data  bases  would  be  needed  for  the  new 
independent  variables  (i.e.  manpower,  weapon  system). 

Further  Research.  As  with  most  research,  a  final  and 
definitive  answer  is  rarely  reached,  and  this  thesis  proved 
*  no  different. 
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The  research  conducted  in  this  thesis  effort  only 
"opened  the  door"  to  finding  a  significantly  accurate  and 
valid  model  for  predicting  SDT  tonnage.  In  statistically 
invalidating  AFLC/DSXR's  current  model,  this  research  has 
caused  an  immediate  need  for  an  accurate  and  valid  model. 
Further  research  in  this  area,  particularly  in  econometric 
models  as  discussed  above,  would  prove  beneficial.  By 
employing  a  model  with  other  independent  variables  such  as 
manpower,  weapon  system  type,  etc.,  as  well  as  variables  tc 
predict  the  unexpected  increases  r  decreases  in  tonnage 
caused  by  diversions,  a  significai  ..y  more  accurate  model 
could  be  developed.  Another  advantage  of  econometric  models 
is  the  ability  to  engage  in  simulation  of  potential  changes 
to  assess  their  impact  on  SDT  budget  requirements. 

Another  aspect  that  may  be  worth  pursuing  is  to  use 
monthly  data  instead  of  quarterly  data.  In  utilizing 
monthly  data,  AFLC/DSXR  will  have  a  far  greater  data  base  to 
work  with  than  the  one  it  currently  uses.  This  could 
improve  the  accuracy  of  any  forecast. 

Further  research  should  also  be  conducted  concerning 
the  idea  to  further  divide  the  SDT  tonnage  being  forecasted 
into  more  definable  subcategories  such  as  particular  weapon 
system  spare  parts  (i.e.  F-16  spares).  This  sub  -it  goriza- 
tion,  although  causing  a  greater  number  of  forecasts  to  be 
conducted,  may  allow  for  more  accurate  estimates  of  SDT 
tcrmaae  . 


88 


Data 


Appendix  A:  SDT  Tonnage  and  Flying  Hour 
SDT  Tonnage  Shipped  to  USAFE  by  MAC 


Quarter 

SDT  Tons 

Flyirg  Hours 

78/1 

9,828 

57 , 480 

/  2 

8 , 942 

52,034 

/  3 

9,793 

64 , 399 

/  4 

9,430 

69 , 807 

79/1 

8  ,  821 

55,423 

/  2 

8 ,831 

56 , 900 

/  3 

8 , 841 

71 , 221 

/  4 

9 , 322 

71 , 846 

80/1 

8 , 496 

60 , 461 

/2 

8 , 573 

58 , 808 

/  3 

8  ,  867 

69 ,735 

/  4 

9 , 174 

68 , 069 

81/1 

9,244 

63 , 546 

/  2 

8  ,  372 

65,588 

/  3 

9,170 

75,263 

/  4 

8 , 747 

75 , 477 

82/1 

8 , 610 

61,878 

/  2 

8 , 138 

70 ,025 

/  3 

9,270 

80 , 973 

/  4 

9,153 

78,365 

83/1 

8,637 

66,217 

/  2 

8 . 242 

70,496 

/  3 

9,178 

77 , 627 

/  4 

9,327 

77 , 944 

84/1 

8,704 

69,485 

/  2 

9,249 

73,093 

/  3 

10 , 058 

78,651 

/  4 

9 , 672 

77 ,702 

85/1 

9 , 289 

69,027 

/  2 

9 , 439 

70 , 119 

/  3 

10 , 259 

83 , 336 

/  4 

9  ,  377 

79,858 

86/1 

9,887 

75 , 552 

/  2 

9 , 383 

74 , 412 

/  3 

10, 508 

80,673 

/  4 

10,700 

77 , 628 

87/1 

8 , 551 

75 , 308 

/  2 

8 , 538 

71 , 368 

/  3 

10, 402 

87 , 519 

/  4 

9 , 602 

81 , 329 

88/1 

9,774 

73 , 531 

/  2 

8 , 254 

71 , 959 

/  3 

6 , 996 

79 , 331 

/  4 

6  ,  848 

83 , 366 

89/1 

7,533 

77 , 571 

/  2 

5,765 

77 , 597 
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SDT  Tonnage  Shipped  to  USAFE  by  MSC 


Quarter 

SDT  Tons 

Flying  Hours 

78/1 

46 , 778 

57 , 480 

/  2 

38 ,450 

52,034 

/  3 

47 , 078 

64 , 399 

/  4 

37 , 450 

69,807 

79/1 

48 , 574 

55,423 

/  2 

44 , 156 

56 , 900 

/  3 

44 , 041 

71 , 221 

/  4 

44 , 296 

71,846 

8C/1 

42 , 971 

60 , 461 

/  2 

50 , 901 

58 , 808 

/  3 

56 , 271 

69,735 

/  4 

49 , 244 

68 , 069 

81/1 

53, 377 

63,546 

/  2 

53,945 

65,588 

/  3 

60 ,785 

75,263 

/  4 

54,941 

75 , 477 

82/1 

55,855 

61 , 878 

/2 

57 , 543 

70,025 

/  3 

63,076 

80,973 

/  4 

65,851 

78 , 365 

83/1 

91 ,436 

66,217 

/  2 

93,263 

70,496 

/  3 

83,737 

77 , 627 

/  4 

83 , 617 

77 , 944 

84/1 

88,315 

69,485 

/  2 

86,968 

73 , 093 

/  3 

101 , 701 

78 , 651 

/  4 

85,521 

77 ,702 

85/1 

96 , 200 

69 , 027 

/  2 

71 , 083 

70 , 119 

/  3 

83 , 702 

83,336 

/  4 

77 , 001 

79,858 

86/1 

75,830 

75,552 

/  2 

79,563 

74 , 412 

/  3 

71,583 

80 , 673 

/  4 

55 , 248 

77 , 628 

87/1 

57 , 088 

75,308 

/  2 

63 , 014 

71 , 368 

/  3 

68 , 675 

87 , 519 

/  4 

69 ,487 

81 , 329 

88/1 

70 , 569 

73 , 531 

/  2 

70,459 

71 , 959 

/  3 

81 , 479 

79 , 331 

/  4 

76 , 619 

83 , 366 

89/1 

73 , 511 

77 , 571 

/  2 

62 , 455 

77 , 597 
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SDT  Tonnage  Shipped  to  PACAF  by  MAC 


Quarter 

SDT  Tons 

Flying  Hours 

78/1 

5,421 

35 , 287 

/  2 

5,801 

37 , 575 

/  3 

6,021 

34,865 

/  4 

5,754 

35 , 922 

79/1 

5,427 

35 , 310 

/  2 

5,456 

37, 593 

/  3 

5,692 

35,600 

/  4 

5,940 

35,615 

80/1 

6 , 495 

35,622 

/  2 

6,647 

37 , 403 

/  3 

6 , 351 

36 , 133 

/  4 

6 , 406 

34 , 429 

81/1 

5,902 

36,236 

/  2 

6,016 

35 , 292 

/  3 

6 , 166 

36 , 480 

/  4 

6,818 

36 , 713 

82/1 

6,363 

36,194 

/  2 

5 , 860 

37 , 007 

/  3 

5,934 

38 , 635 

/  4 

5,368 

37 , 534 

83/1 

5,729 

37 , 293 

/  2 

6,768 

39,678 

/  3 

6 , 386 

39,129 

/  4 

6,203 

37 , 617 

84/1 

6 , 020 

40 , 018 

/  2 

6 , 916 

40 , 533 

/  3 

6 ,050 

39,523 

/  4 

5,829 

38 , 235 

85/1 

5,792 

40 , 802 

/  2 

6 ,188 

40 , 828 

/  3 

5,938 

41,344 

/  4 

6,084 

40 , 983 

86/1 

5,829 

42 , 905 

/  2 

5,569 

41 , 942 

/  3 

6,782 

41 , 476 

/  4 

6 , 285 

41,372 

87/1 

5 ,830 

44 , 270 

12 

6 , 008 

42,322 

/  3 

6  ,  566 

44 , 381 

/  4 

5,886 

43 , 700 

88/1 

5,835 

42,588 

/  2 

5  ,  582 

41 , 962 

/  3 

4,787 

42 , 206 

/  4 

3,906 

42,255 

89/1 

3,850 

42,856 

/  2 

4 ,039 

43 , 482 
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SDT  Tonnage  Shipped  to  PACAF  by  MSC 


Quarter  SDT  Tons 


Flying  Hours 


78/1 

31,163 

/  2 

30,967 

/  3 

32 ,924 

/  4 

32 , 018 

79/1 

33 , 469 

/  2 

30,548 

/  3 

33,046 

/  4 

34 , 991 

80/1 

33,145 

/  2 

30 , 312 

/  3 

35 , 918 

/  4 

32,220 

81/1 

35,198 

/  2 

30 , 649 

/  3 

35,193 

/  4 

35,396 

82/1 

37 , 343 

/  2 

41,379 

/  3 

43,392 

/  4 

42,968 

83/1 

43 , 039 

/  2 

49,651 

/  3 

46 , 352 

/  4 

35,398 

84/1 

38 , 462 

/  2 

41 , 800 

/  3 

48 , 352 

/  4 

49 , 203 

85/1 

47,567 

/2 

49, 835 

/  3 

59, 435 

/  4 

48 , 235 

86/1 

49 , 398 

/  2 

41,209 

/  3 

41,782 

/  4 

34 , 195 

87/1 

41 , 814 

/  2 

39 , 046 

/  3 

34,185 

/  4 

36,701 

88/1 

42 , 387 

/  2 

48 , 394 

/  3 

55,113 

/  4 

42 , 250 

89/1 

42 , 513 

/  2 

44 , 280 

35,287 
37 , 575 
34 , 865 
35,922 
35, 310 
37 , 593 
35,600 
35,615 
35,622 
37 , 403 
36 , 133 
34,429 
36,236 

35.292 
36 , 480 
36 ,718 
36 , 194 
37 , 007 
38,635 
37 , 534 

37 . 293 
39,678 
39,129 
37,617 
40 ,018 
40 , 533 
39,523 
38 ,235 
40 , 802 
40,828 
41,344 
40 , 983 
42 , 905 
41,942 
41 , 476 
41,372 
44 , 270 
42 , 322 
44 , 381 
43  ''OO 
42, 5?S 
41,962 
42 , 206 
42 , 255 
42,856 
43,482 
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Appendix  B 


Flying  Hour  Parameter  Confid ence 
Interval  Values 


SDT 

Tonnage  Shipped 

to  USAFE  by 

MAC 

#  of 

Lower 

Upper 

Quarters 

Limit 

Parameter 

Limit 

40 

0.0173 

0 .0385 

0 .0598 

39 

0.0241 

0 .0449 

0.0^58 

38 

0 .0283 

0.0507 

0.0732 

37 

0.0326 

0.0544 

0.0763 

36 

0 .0329 

0.0549 

0 . 0769 

35 

0 .0368 

0 . 0604 

0 .0839 

34 

0 . 0416 

0 . 0668 

0 . 0920 

33 

0 .0410 

0 . 0665 

0.0920 

32 

0.0405 

0.0667 

0. 0929 

31 

0.0398 

0.0681 

0.0964 

30 

0 . 0417 

0 . 0731 

0 . 1044 

29 

0.0405 

0.0728 

0 . 1051 

28 

0.0417 

0 . 0751 

0 . 1084 

27 

0.0508 

0.0852 

0.1195 

26 

0.0469 

0.0839 

0.1210 

25 

0.0461 

0.0839 

0.1218 

24 

0.0463 

0.0840 

0.1217 

23 

0 . 0497 

0 . 0931 

0 .1365 

22 

0.0421 

0 .0859 

0 . 1297 

21 

0.0465 

0.0909 

0.1354 

20 

0.0484 

0.0932 

0.1380 

19 

0.0^41 

0.0951 

0.1460 

18 

0 .0333 

0.0849 

0 .1364 

17 

0.0339 

0 .0859 

0 .1380 

16 

0 .0339 

0 .0872 

0 . 1404 

15 

0 . 0225 

0 . 0814 

0.1402 

14 

0.0171 

0.0803 

0 . 1434 

13 

0.0131 

0 . 0792 

0 .1453 

12 

0 .0090 

0 . 0792 

0 .1494 

11 

0 . 0049 

0 . 0881 

0 . 1713 

10 

0 . 0068 

0 .1061 

0 . 2053 

9 

-0.0109 

0.1047 

0 . 2204 

8 

-0 .0156 

0.1096 

0 .2347 
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SDT  Tonnage  Shipped  to  USAFE  by  MSC 


#  of 

Quarters 

40 

39 

38 

37 

36 

35 

34 

33 

32 

31 

30 

29 

28 

27 

26 

25 

24 

23 

22 

21 

20 

19 

18 

17 

16 

15 

14 

13 

12 

11 

10 

9 

8 


Lower 

Limit 

Parameter 

Upper 

Limit 

0 .5329 

1 . 1275 

1 . 7220 

0.4869 

1.1117 

1 .7365 

0 . 3729 

1 . 0607 

1 .7484 

0.3165 

1 . 0188 

1 .7210 

0.3045 

0 . 9868 

1 . 6692 

0.2165 

0 . 9641 

1.7116 

0 .0530 

0.8703 

1 . 6876 

0 . 0370 

0.8390 

1.6411 

0 . 0257 

0.8157 

1.6057 

-0.1867 

0 .6469 

1 . 4805 

-0 .4271 

0 .4973 

1 . 4216 

-0 .4982 

0.4397 

1 . 3775 

-0 . 6272 

0 . 3146 

1.2565 

-0 . 8906 

0 .1078 

1.1063 

-1.1546 

-0 .1202 

0.9141 

-1 .1592 

-0 . 1219 

0.9153 

-1.1324 

-0 . 1190 

0 .8945 

-1.7070 

-0.6126 

0 .4819 

-1 . 9169 

-0 .8562 

0.2046 

-1.8531 

-0.7614 

0 . 3303 

-1 . 8197 

-0.7125 

0.3947 

-1 . 8447 

-0. 5923 

0.6602 

-1 .7416 

-0 . 4259 

0 .8898 

-1.7966 

-0.4409 

0 . 9147 

-1 . 8650 

-0 . 4636 

0 .9379 

-1 .8387 

-0.2957 

1.2473 

-1 . 8011 

-0 . 1787 

1 . 4437 

-1 . 6941 

-0 . 2887 

1 . 1167 

-1 .7191 

-0.3073 

1 . 1045 

-C  .  8706 

0 . 3870 

1 . 6447 

-0 . 9931 

0 . 5492 

2 . 0914 

-1 . 3638 

0 . 2371 

1  .  8380 

-1 . 5351 

0 .1530 

1 . 8411 
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SDT  Tonnage  Shipped  to  PACAF  by  MAC 


#  of 

Lower 

Quarters 

Limit 

40 

-0 .0360 

39 

-0 . 0430 

38 

-0 .0445 

37 

-0.0473 

36 

-0 .0520 

35 

-0 . 0613 

34 

-0.0638 

33 

-0 . 0716 

32 

-0 . 0784 

31 

-0 . 0768 

30 

-0.0732 

29 

-0 .0597 

28 

-0.0567 

27 

-0 . 0641 

26 

-0.0734 

25 

-0 . 0776 

24 

-0 .0576 

23 

-0 .0502 

22 

-0 . 0622 

21 

-0 . 0692 

20 

-0 .0988 

19 

-0 . 1326 

18 

-0 . 1163 

17 

-0. 1137 

16 

-0 . 1290 

15 

-0.1418 

14 

-0 .0984 

13 

-0 .1107 

12 

-0. 1902 

11 

-0 . 2322 

10 

-0 .2597 

9 

-0 . 3078 

3 

-0. 3912 

Parameter 

Upper 
Limi  t 

0 . 0126 

0 . 0612 

0 . 0058 

0 . 0545 

0 . 0047 

0 .0539 

0 . 0038 

0 .0550 

0 . 0002 

0 .0524 

-0 . 0093 

0 . 042S 

-0.0131 

0 .0377 

-0 .0202 

0.0312 

-0 . 0247 

0 . 0290 

-0.0209 

0 . 0351 

-0.0168 

0 .0395 

-0 . 0047 

0 .0503 

0 . 0029 

0 . 0625 

-0 . 0013 

0.0614 

-0 . 0048 

0.0638 

-0 . 0040 

0 . 0695 

0 .0160 

0 .0895 

0 . 0297 

0 . 1097 

0 . 0252 

0 .1126 

0 . 0227 

0 . 1147 

-0 . 0059 

0 . 0870 

-  0 .0313 

0.0700 

-0 . 0174 

0 . 0S15 

-0 . 0086 

0 . 0965 

-0 . 0024 

0 . 1242 

-0 . 0057 

0.1303 

0.0189 

0 . 1362 

0 . 0222 

0 .1550 

-0 . 0058 

0.1787 

-0.0307 

0.1708 

-0 . 0245 

0 . 2106 

-0.0404 

0.2269 

-0 .0536 

0.2340 
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SDT  Tonnage  Shipped  to  PACAF  by  MSG 


#  of 

Quarters 

40 
39 
38 
37 
36 
35 
34 
33 
32 
31 
30 
29 
28 
27 
26 
25 
24 
23 
22 
21 
20 
19 
'  13 
17 
16 
15 
14 
13 
12 
11 
10 
9 
8 


Lower 

Limit 

0 . 6074 
0 . 5  ‘74 
0 . 5230 
0.4685 
0 . 4061 
0 . 3447 
0 . 3061 
0 . 2273 
0 .1535 
0 .0504 
-0 .0248 
-0 .1079 
-0 . 33 37 
-0.4777 
-0 .8120 
-1 .0399 
-1 . 2°  3  3 
-1 .6539 
-1 . 9044 
-2 . 0104 
-2 . 2625 
-2.6618 
-2.6172 
-2 .7108 
-3 . 8217 
-4 .1594 
-4 . 3907 
-4 .7401 
-6 . 3182 
-6 . 7867 
-7 .1545 
-4 . 9062 
-4 .4304 


parameter 

1 . 3024 
1 . 2609 
1 .2346 
1 . 2078 
1 . 1621 
1 .1290 
1 .0848 
1 .0326 
1 . 0008 
0.9324 
0 . 8499 
0 .8046 
0.6494 
0 . 5522 
0 . 261 5 
0.0812 
-0.1342 
-0 . 4025 
-0.5466 
-0 . 5812 
-0 . 6968 
-0 .8925 
-0 . 7809 
-0.7333 
-1.7134 
-2 . 0360 
-2 . 1780 
-2 . 2334 
-2 .8553 
-2.8574 
-2 . 5967 
-1 . 2747 
-0 . 3093 


Upper 

Limit 

1 . 9975 
x  .  9744 
1 .9461 
1 .9470 
1.9180 
1 . 9133 
1 .8636 
1 .8378 
1 .8482 
1.8144 
1 .7245 
1 .7172 
1 . 6325 
1 . 5821 
1 .3351 
1 . 2023 
1  0299 
0 . 8486 
0.8113 
0 . 8481 
0 . 8690 
0 . 8767 
1 . 0554 
1 . 2442 
0 . 3949 
0 . 0875 
0.0343 
0.2732 
0 . 6076 
1.0719 
1 . 9610 
2 .3568 
3 .8117 
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Appendix  C:  Original  Series  Autocorrelation  and  Partial 
Autocorrelation  Functions 


AUTOCORRELATION  FUNCTION 

DATA  -  SDT  Tonnage  Shipped  ay  SAC  ::  'JS Arc 

DIFFERENCING  -  ORIGINAL  SERIES  IS  YOUR  DATA. 

DIFFERENCES  5SL0M  ARE  OF  ORDER 


GRIGISAL  SERIES 

KEAN  OF  THE  SERIES  =  .?0I56E*D4 
ST  DEV  rr  ‘JEP7”  * 

SUSSES  Or’OBSEEVATIONS”  ”“40 


24 


.24  .24 


.05 

’  i 


SEAN  DIVIDED  3:  ST.  ERROR  -  ,53162E»02 

TO  TEST  iHETSER  THIS  SERIES  IS  WHITE  NOISE.  THE  VALUE 
SHOULD  3E  COMPARED  WITH  A  OHI -SQUARE  VARIABLE  WITH  2 


■  c 

.  L  I 


■  i  i 

.  1  L 


.  4i9:2E*02 
DEGREES  CF  r 


.  m  1  .  i  1 


DIFFERENCE  I 

MEAN  0?  THE  SERIES  =- . 5T949E-0 


cfsrrc  : 


.  .  1  -Si*.' 


NUMBER  CF  OBSERVATIONS  - 

::  12  -.11  -.43  .03  -.03  -.26  .00 

II-  24  -.14  -.21 

SEAN  DIVIDED  BY  ST.  ERROR.  *  .4S443S-QI 

TO  TEST  WHETHER.  THIS  SERIES  IS  WHITE  NOISE.  THE  VALUE  . 441  j'E* 
SHOULD  BE  COMPARED  WITH  A  CHI-SQUARE  VARIABLE  WITH  24  DEGREES  0 


.03 

.26 


SEAN  Or  THE  SERIFS  =  .12632 

j  *  <J  1 

’’l  p  ;  *  '  J»  »  >  •  •»  g  .  4 

mivcra  r. rct5 = 

3d 

:  .  i "  t  « 

*  *  id  .  t : 

.2"  -.33 

v 

s”  ;  ’ s  * : 

t  r. 

Ar.  y. 

*  * 

2’  y 

13*  24  -.12  -.15 

~ ,  v  4 

.36  -.22 

-.03 

_  ■u  A ;  .  i 

.  v  J  •  4  -  .  *  V 

3T.E.  .24  .24 

.24 

.24  .16 

.  i! 

.26  .26  .2' 

MEAN  DIVIDED  BY  ST.  ERROR  = 

.  ;i5j 

9--G1 

mQ  I’rc’1  tdt<  crp’jc 

IS  WHI 

TE  NOISE,  TH 

E  VALUE 

■  •BeLOE-DI 

SHOULD  BE  COMPARED  WITH  A  I 

HI-SQ-JA 

RE  VARIABLE 

WITH  24 

DEGREES  OF  FREEDOM 
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SDT  Tannace  Shipped  by  MAC  cc  VSAFE 
GRAPH  CP  Dl FFERtMCE  2  ACF 

5.000  10.000  15.000 


M  ,inr 


25.000 


900 


750 


500 


Hi 


300 


’50 


.000 


2  50  + 


J'JO 


/* 

0  0  0  * 

0  X  : 

0  *  0 

0  0 

5  * 

2  5 

0  0  * 

i,,1 

j  2  •  * 

X 

X  X 

X 

X  X 

X 

7  7 

X 

X  X 

X 

X  “  * 

X 

y 

7 

X  X 

y 

7 

X 

X 

1 

X 

V 

X 

X 

X  X 

X 

Y 

7 

v; 

X 

7 

Y 

Y 

V 

v 

X 

X 

Y 

X 

X 

X  + 

X 

+ 

X  * 

X 

* 

X 

V 

X 

7 

X 

X 

Y 

7 

x 

X 

X 

X 

X 

X 

X 

X 

Y 

t 

* 

4. 

t 

-Hv 
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PARTIAL  AUTOCCRRELATICNS 

DATA  -  SDT  Tonnage  Shipped  by  MAC  cc  USAFE 

DIFFERENCING  -  ORIGINAL  SERIES  IS  YOUR  DATA. 

DIFFERENCES  BELOW  ARE  Of  ORDER 


ORIGINAL  SERIES 

MEAN  OF  THE  SERIES  =  .92156S+G4 
ST.  DEV.  OF  SERIES  =  . 5D56jE*G5 
N'JMBER  OF  OBSERVATIONS  =  40 


DIFFERENCE  I 
MEAN  OF  THE 


NUMBER  IF  OBSERVATIONS 


ST.  DEV.  Cf  SE5!!:  :  .III 
NUMBER  CFVOBSERVATIONS“’ 


99 


*** "  3T  a  *  °  SnivC"d  CV  v" 

■“-L'j  efi'JF^Isq^Ficf 


;:  USA: E 


*50X0 
X  X 
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AUTOCORRELATION  FUNCTION 

DATA  -  SDT  Tonnage  Snipped  by  KSC  to  U SAFE 

uiFFERENCi.rU  -  ORISL.Au  SERIES  IS  YOUR  DATA. 

DIFFERENCES  BELOV  ARE  OF  ORDER 


ORIGINAL  SERIES 

KEAN  OF  THE  SERIES  =  .54980E+05 
ST.  DEV.  OF  SERIES  =  . 17731E+D5 
NUMBER  OF  OBSERVATIONS  =  40 


FRIES  IS  WHITE  NOISE.  THE  VALUE  .I'-ISE-OJ 


SHOULD  BE  COMPARED  KITH  A  CHI-SQUARE  VARIABLE  WITH 


DIFFERENCE  1 

MEAN  OF  THE  SERIES  =  .532232+03 
ST.  DEV.  OF  SERIES  =  .92546E-04 
NUMBER  OF  OBSERVATIONS  =  35 


ii  ucar.;:; 


MEAN  DIVIDED  BY  ST.  ERROR  =  .35252E+03 

TO  TEST  WHETHER  THIS  SERIES  IS  WHITE  NOISE.  THE  VALUE  .23 
SHOULD  BE  COMPARED  WITH  A  CHI-SQUARE  VARIABLE  WITH  24  DEG? 


DIFFERENCE  2 

MEAN  OF' THE  SERIES  =  .24053E+03 
ST.  DEV.  OF  SERIES  =  . I5534E+05 
NUMBER  OF  OBSERVATIONS  =  33 


MEAN  DIVIDED  BY  ST.  ERROR  =  .55443E-01 

TO  TEST  WHETHER  THIS  SERIES  IS  WHITE  NOISE,  THE  VALUE  .56162E+02 
SHOULD  BE  COMPARED  WITH  A  CHI-SQUARE  VARIABLE  WITH  24  DEGREES  Or  FREEDOM 


42  05SZRY- 
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.15  .63 

.  5  ‘ 

.56 

. L .  ... 

■ 7 

* 

ST.E.  .16 

n  r  q  ft 
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.35  .‘40 

1  1  »  A 

.  1  .  .N; 

'42 

1  - 
.  1  . 

.  i  - 

3-  24  -.15 

-.23  -.31 

-.31 

—  .  r 

-.33  -.42 

-.44  -.41 

. 

-.56 

ST.E.  .42 

.42  .42 

.43 

A] 

.44  .45 

i  *  i  ** 

.  *2  .t 

•  1  t 

.  1 1 

.  : 

AN  DIVIDED  BY  ST. 

ERROR  =  .231 

7  5  E + 0  2 
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PARTIAL  AUTOCORRELATIONS 

DATA  -  SDT  Tonnage  Shipped  by  NSC  to  USAFE 

DIFFERENCING  -  ORIGINAL  SERIES  IS  YOUR  DATA. 

DIFFERENCES  BELOW  ARE  OF  ORDER  1 


ORIGINAL  SERIES 

MEAN  OF  THE  SERIES  =  .54980E+Q5 
ST.  DEV.  OF  SERIES  =  .17731E+05 
NUMBER  OF  OBSERVATIONS  --  40 

I-  12  .35  .24  -.15 

13-  24  -.05  -.91  -.13 


DIFFERENCE  1 


MEAN  OF  TEE 

SERIES  = 

.53223**93 

ST.  DEV.  OF 

SERIES  = 

.  3254SE+34 

NUMBER  OF  OBSERVATION 

)S  =  33 

i.  n 

•  .  u 

-.33 

.12  -.15  -.39 

.04  .21  .3? 

13-  24 

-.03 

.05  -.24  .0? 

.03  -.23  -.31 

DIFFER ^NCE  5 

MEAN  OF- THEORIES  =  .240535*33 
ST.  DEV.  OF  SERIES  =  . 15534E+05 
NUMBER  OF  03SERVATICNS  =  33 

1-  12  -.72  -.24  -.24 

13-  24  -.03  .14  -.15 
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SDT  Tonnaae  Shipped  hy  MSC  t:  USAFE 
GRAPH  OF  DIFFERENCE  1  FACF 
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AUTOCORRELATION  FUNCTION 

DATA  -  SDT  Tonnage  Shipped  by  MAC  :c  PACAF 

DIFFERENCING  -  ORIGINAL  SERIES  IS  YCJR  DATA. 

DIFFERENCES  BELOW  ARE  OF  ORDER  1 


ORIGINAL  SERIES 

MEAN  OF  THE  SERIES  =  .60774E+04 
ST.  DEV.  OF  SERIES  =  . 42359E+03 
NUMBER  OF  OBSERVATIONS  =  40 
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.09 
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13-  24 

.23 
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A  *> 
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•  U  J  .  U  v 

.20 

.20 

!H  .’n 

A  f 
.  L  1 

MEAN  DIVIC 

ED  BY  ST. 

ERROR  = 

.90742E+02 

TO  TEST  WH 

ETHER  THIS 

SERIES 

IS  WHITE  NCISE 

THE 

VALUE 

.153535+02 

*  w  i  U  J  i  +  Jwl\*uw  it,  n.liiu  k.i*  lUDk  *  u 

SHOULD  BE  COMPARED  WITH  A  CHI-SQUARE  VARIABLE  KITH  24  DEGREES  OF  FREEDOM 


DIFFERENCE  1 

MEAN  OF  THE  SERIES  =  .119235+32 
ST.  DEV.  OF  SERIES  =  .4S217E+03 
NUMBER  OF  OBSERVATIONS  =  j? 
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MEAN  DIVIDED  BY  ST.  ERROR  --  .151435+90 


TO  TEST  WHETHER  THIS  SERIES  IS  WHITE  NOISE.  THE  VALUE  .159S7E+92 
SHOULD  BE  COMPARED  WITH  A  CHI-SQUARE  VARIABLE  WITH  24  DEGREES  OF  FREEDOM 


DIFFERENCE  2 

MEAN  OF  THE  SERIES  -  .273355+02 
ST.  DEV.  OF  SERIES  =  .739005+03 
NUMBER  OF  OBSERVATIONS  =  33 
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MEAN  DIVIDED  BY  ST.  ERROR  =  . 23253E+00 

TO  TEST  WHETHER  THIS  SERIES  IS  WHITE  NOISE.  THE  VALUE  . 19272E+32 
SHOULD  BE  COMPARED  KITH  A  CHI-SQUARE  VARIABLE  WITH  24  DEGREES  OF  FREEDOM 
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SDT  Tcnnaas  Shipped  by  MAC  to  PACAF 
GRAPH  OF  OBSERVES  SERIES  ACf 
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PARTIAL  AUTOCORRELATIONS 

DATA  -  SDT  Tonnage  Shipped  by  MAC  :o  PACAF 

DIFFERENCING  -  ORIGINAL  SERIES  IS  YOUR  DATA. 


DIFF 

ERENCES  BELOV 

ARE  OF 

ORDER  2 

ORIGINAL  SERIES 

MEAN  OF  THE 

SERIES 

=  .5077  4£»04 

ST.  DEV.  OF 

SERIES 

=  .42353E+93 

NUMBER  OF  OBSERVATIONS  =  40 

I  12 
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'}0 

-.22  .05 

13-  24 
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-.13  -.14 
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NUMBER  OF  0 
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1 1 
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ST  .  DEV.  OF 
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SOT  Tcnnacs  Shipped  by  MAC  tc  PACAF 

graph  c?  Observed  series  pace 
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AUTOCORRELATION  FUNCTION 

DATA  -  SDT  Tonnage  Chipped  by  NSC  to  PACAF 

DIFFERENCING  -  ORIGINAL  SERIES  IS  YOUR  DATA. 

DIFFERENCES  BELOW  ARE  OF  ARDER  1 

ORIGINAL  SERIES 

MEAN  OF  THE  SERIES  =  . 33193E+05 
ST.  DEV.  OF  SE?TES  -  .?C'50£*J4 
NUMBER  OF  OBSERVATIONS  =  40 


DIFFERENCE  C 


m  £  a  N  OF  THE  SERIES  r  . “13:5: ♦02 
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Ill 


3D7  Tcnnaie  Shipcei  fcv  “00  +:  PACAr 

graph  of  Difference  i'paof 
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Appendix  D:  TIMES  Forecasting  Computer  Output 

SUMMARY  OF  5DT  Tonnage  Shipped  by  MAC  tc  ’'SAFE 

(tmttninitunittttttittitittttttttt  ttttntt  tHMtiimntunimitnimtittittmtmtiittnttii 

DATA  -  2  =  SDT  Tonnage  Shipped  by  MAC  to  USAFE  13  OBSERVATIONS 

DIFFERENCING  ON  2  -  2  OF  ORDER  1 

1 1 1  i  1 1 1  (  n  1 1  n  1 1 1 1 1 1  ii  1 1 1 1  n  u  m  1 1 1 1 1  1 1  u  1 1 1  ■  i  n  n  t  n  1 1  n  i  i  u  u  i  n  m  t  1 1 1  i  n  m  i  1 1  i  i  1 1 1  i  1 1  r  r  1 1  m  m  u 


PARAMETER  PARAMETER 

NUMBER  TYPE 


PARAMETER  r  S  T I  MATED 

"  order"'  ’  value'" 


o;  :tb  -:■>;* 

/  J  4  .Ull. 


LOVER  LIMIT 

ummiutitttktutmmumttttimntuuiiitntuunutututtui  tmttiuim  t  iiitttiituittu 

1  AUTOREGRESSIVE  I  2  -.oC-TISnO  -.31116EOO  - . ?"E-L: 

2  MOVING  AVERAGE  1  1  .IITRIE^I  .  1  14552*::  . I'll EE*:i 

3  MOVING  AVERAGE  I  2  -.T31T4S+00  -.I0229E*01  - . 4 4 G 5 5 E - 0 1 

tttttmitntttuittitmntttumnttttuutHtttuitmmuttitttititimtttttittittiittiiittttiiii 

OTHER  INFORMATION  AND  RESULTS 

ittitiittittttttttiuituittmitttmtimuuttttt  ununtuutHt  ttu  tuuttttut  it  tiiuit  tit  in  it  t 

RESIDUAL  SUM  Or  SQUARES  .U133E*0S  33  D.F .  RESIDUAL  MEAN  SQUARE  .3:';2E-:6 

NUMBER  OF  RESIDUALS  36  RESIDUAL  STANDARD  ERROR  .6323:1-:: 
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THE  ESTIMATED  RESIDUALS  --  MODEL  I 
OBSERVED  SERIES 
DEVIATIONS'"'??.?!!  THE  MEAN 

11.300  21.000  31.000  41.000  51.000  61.000  ’1.000  31.000  91.000 
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AUTOCORRELATION  FUNCTION 

DATA  -  THE  ESTIMATED  RESIDUALS  --  SDT  Tonnage  Shipped  by  MAC  :: 
DIFFERENCING  -  ORIGINAL  SERIES  ;S  '/OUR  DATA. 

DIFFERENCES  BELOW  ARE  Of  ORDER  I 

ORIGINAL  SERIES 

MEAN  OF  THE  SERIES  =  J7853E+02 
ST.  DEV.  OF  SERIES  =  .61D54E+Q3 
NUMBER  OF  OBSERVATIONS  =  3: 
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MEAN  DIVIDED  3/  ST.  ERROR  =  ,47024E-‘00 

"lfl  TrST  WHETHER  THIS  SER*"S  rS  ilHIT*'  NOISE.  THE  V ALU 
S  H  O  U  L  D  ‘  3  E  *  ‘  C  D  M  F  A  E  E  D  *  W I T  H  ‘  '.A "  C  H I  -  S  -C  U  A  R  E  VARIABLE  VITh"" 


'.if1 Ji 
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THE  ESTIMATED  RESIDUALS  -  SET  Tannage  Shipped  bv  MAC  :c  USAFE 
GRAPH  CF  OBSERVED  SERIES  ACF 
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THE  ESTIMATED  RESIDUALS  -  SDT  T:nnaae  Shipped  by  MAC  to  L'SAFE 
CUMULATIVE  PERIODOCRAM  .1  PROBABILITY  LIMITS 
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TEE  ESTIMATED  RESIDUALS  -  SDT  Tonnage  Shipped  bp  MAC  to  L'SAFS 
HISTOGRAM 
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PREKHITEIE3  SDT  Tcsnace  Shipped  by  MAC  to  USAPE 

LQGl 0  SPECTRUM  SMOOTHING  3AND¥ luTH  =  .393  APPROX  35  P.C.  CONFIDENCE  LIMITS 
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DATA  -  l  -  SDT  Tonnage  Shipped  by  MAC  to  USAFE  10  OBSERVATIONS 

DIFFERENCING  ON  2  -  2  OF  ORDER  1 

PARAMETER  PARAMETER  PARAMETER  ESTIMATED 

NUMBER  TYPE  ORDER  VALUE 

tlilliilimimlimtiitiiiiiittitiiiiiiiiiiitiimittiiiiiiiitiiitiiiiittumititiiimiiitimiiin 

1  AUTOREGRESSIVE  I  2  -.50572E+00 

2  MOVING  AVERAGE  1  1  .14292E+01 

3  MOVING  AVERAGE  1  2  -.73I74E+00 

mimiltllittilttlittlnliititmintltlmtiitmiilimitmilinimtinimlmiiitiitiiiiiimti 

NUMBER  OF  TIME  ORIGINS  FOR  FORECASTS  -  1 

NUMBER  OF  FORECASTS  AT  EACH  TIME  ORIGIN  =  6 

FORECAST  TIMS  ORIGINS  ARE  T  =  4Q 

SET  Tonnage  Shipped  by  MAC  to  USAFE  FORECASTS  BASE  PERIOD  40  WITH  95  PER  CENT  CONFIDENCE  LIMITS 
PERIODS  AHEAD  LO.  CONF.  LIMIT  FORECAST  UP.  CONF.  LIMIT 

1  .7279630E*04  .85IS104E*34  .9752579E^C4 

2  .7207063E+04  .3630789E+34  . 10054515*05 

3  .  7  395529E+C4  .88900575*04  .  1038458E*05 

4  .6683235E+04  .8542168E+04  .10396104+05 

5  . 5673475E+04  .8121161E+34  .10568S4E*05 

6  . 5056806E+04  .30072055+04  .lu95760E*05 


SUMMARY  OP  SDT  Tonnage  Shipped  by  MSC  t:  USAFE 

iittiiiuniittiuitttiittuiitutiuttttiitutiiuttuttttiutiimtnitinmttmnttittmitmimt 

DATA  -  S  =  SDT  Tonnage  Shipped  by  MSC  to  USAFE  43  G3SER7ATIC.VS 

DIFFERENCING  CN  2  -  1  OF  ORDER  1 

ittmitituttttttttttttttumtitiuumuuttHUiiitttttuittitttitmttttiinuttittuitntttt  tut 

PARAMETER  PARAMETER  PARAMETER  ESTIMATED  95  PER  CENT 

NUMBER  TYPE  ORDER  VALUE  LONER  LIMIT  UPPER  LIMIT 

itttitttttitttntttttnitiittuitimttntntiiitutnimtmtntmtuun:  tittuttutmiuii  ttt  mi 


1 

AUTOREGRESSIVE  1 

1 

- . SC9OQE+0Q 

-.11324E+01 

- .  4  3  5  6 1 E  t  C-  0 

l 

MOVING  AVERAGE  I 

1 

-  5753flE*00 

1I220E-C1 

- . 25526E“*1 

ttttttttAiiinmtnHnniiJtiitiitiiititititiiimtiiittitttittttittitt  tittintn  it  mi  tin  tut  in  tin 


OTHER  INFORMATION  AND  RESULTS 

inttumttuttikmtnimtimttttititMittiimntiuuttumiitittiimtittmtnimtutituuiui 


RESIDUAL  SUM  OF  SQUARES 

;5«' triin 

<  U  W  .  *.  -J  M  i  w 

36  D.F. 

RESIDUAL  MEAN  SQUARE 

."15c"E*C3 

NUMBER  OF  RESIDUALS 

33 

RESIDUAL  STANDARD  ERROR 

- 34345E*C4 
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THE  ESTIMATED  RESIDUALS  -  SDT  Tonnaae  Shipped  by  MSC  to  USA? E 
OBSERVED  SERIES 
DEVIATIONS  FROM  THE  MEAN 

11.009  21.000  31.000  41.000  51.000  61.000  71.900 


81.000  91.300 
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THE  ESTIMATED  RESIDUALS  --  SDT  Tcn33',$  Siirr$d  bv  MSI  U:A:? 
CUMUjATI V£ " PER i5D2GRAM^  .1  ? ?. 0 5 A B  ■  L  ■  f Y^L ! M : ^ *S * r ^ . **’ 
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PARAMETER  PARAMETER  PARAMETER  ESTIMATED 

NUMBER  TYPE  ORDER  VALUE 


AUTOREGRESSIVE  1  1  30900E-00 

MOVING  AVERAGE  1  1  57530E+0C 


NUMBER  0:  TIME  ORIGINS  FOR  FORETASTE  = 
NUMBER  OF  FORECASTS  AT  EACH  TIME  ORIGIN  = 


CAST  i IME  ORIGIN 

S  ARE  T  = 

40 

Tonnage  Shipped 

by  MSC  ::  USAFE 

FORECASTS  BASE  PERIOD 

40  iilTH  55  PER  CENT  OOSFIDENOE  LIMITS 

CDS  AHEAD  LO 

.  CON?.  LIMIT 

FORECAST 

UP.  TON:.  LIMIT 

f 

.5339166E+05 

.'002123E+05 

.  3  5  5  5  C1 9 1 2  *  0  5 

*. 

.4863526E-05 

.  6  9  5  3  9  0  5  E  -  0  5 

Q '  5  2  9  S  4  f  - 

< 

, 4364533E+05 

.59933732*05 

!962  2213E^5 

4 

.401TC25E+05 

.  6  9555S4E^05 

.  99141422-35 

* 

.3551 300E+05 

. 59S84V0S+35 

.10315042-06 

r\ 

.  337C913S+0; 

.59699552-05 

. 1T55B99E-L6 
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smm  OF  507  Teams  Shipped  by  MAC  to  FACAF 


DATA  -  2  =  SDT  Tonnage  Shipped  by  MAC  tc  FACAF  40  C::i?.7A7I0!iS 

DIFFERENCING  ON  l  -  NONE 


THE  ESTIMATED  RESIDUALS  -  SDT  Tennace  Shipped  by  MAC  r:  USAS: 

OBSERVED  SERIES 
DEVIATIONS  FROM  THE  MEAN 

11.000  21.300  31.000  41.000  51.000  61.000  ’’1.000  31.303 


X  X 


X  X  X 


+  -  -  + 


.433  +  -X+ 


.003  *XMMMMMMMM*MMXMMXMMM+MXMMMMMMM*MMNMMMXMM»MKMMMMMMM-'MMMMMMMMM’MMMMMMMMMiMMMMMMMMM*MMMMMMMMM*MMMMMMMMM* 


.230  *  ---X  ‘X  X  -+  -  -  *  --  X  -* 

-  ...  -  -XX  y 


-.430  *  -X  * 


X  *  X  x+ 
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AUTOCORRELATION  FUNCTION 
DATA  -  THE  ESTIMATED  RESIDUALS 


SDT  Tcnr.age  Shipped  by  MAC  FACAF 


DIFFERENCING  -  ORIGINAL  SERIES  IS  YOUR  DATA. 

DIFFERENCES  BELOV  ARE  OF  ORDER  1 


ORIGINAL  SERIES 

MEAN  OF  THE  SERIES  =  .E5003E‘01 
ST.  DEV.  OF  SERIES  =  .3?25:E+jJ 
NUMBER  OF  OBSERVATIONS  =  40 

I-  12  .00  -.04  -.07 

S  tfi  i:  .In 


MEAN  DIVIDED  BY  ST.  ERROR  -  .1B473E*00 

TO  TEST  WHETHER  THIS  SERIES  IS  WHITE  NOISE.  THE  VALUE 
SHOULD  BE  COMPARED  WITH  A  CHI-SQUARE  VARIABLE  WITH  22 


.54246E-C! 

DEGREES  OF  FREEDOM 
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THE  ESTIMATED  RESIDUALS  - 
GRAPH  OF  OBSERVED  SERIES  AC? 

5.303  10.000  15. 


SDT  Tonnage  Shipped  by  MAC  ::  PACAF 
030  23.030  25.000  30.000  35. 


10. 


.050 


.900 


.750 


.600 


.150 


.300 


.150 


nnn 


j  'j 

x  x 

x 


1  SO 


inn  + 


v 

+  0  3 
X  X 
X 
X 
X 


0  0  *  3  0  0 
X  X  X  X  X 
X  X 
X 

V 

X 


X  X 

n 


X 

0  X 


y 

X  X 

n  n  +  n 


+  n 


132 


THE  ESTIMATED  RESIDUALS  -  SDT  Tonasce  Shifted  cv  MAC  cc  PACAF 
CUMULATIVE  PERIODOORAM  .1  PROBABILITY  LIMITS 
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THE  ESTIMATED  RESIDUALS 
HISTOGRAM 

CELL  WIDTH  .500CE+03 
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PREBKITENED  SD?  Tcnnsje  Shipped  bv  MAC  tc  PACAf 

L0510  SFECTHUM  SMOOTHING  BANDWIDTH  =  .CSC  APPROX  35  ?.C.  CONFIDENCE  LIMITS 
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TIME  SERIES  FORECASTING  FOP.  SDT  Tonnage  Shipped  by  MAC  :c  PAG 


DATA  -  2  =  SDT  Tonnage  Shipped  by  MAC  to  PACAF 
DIFFERENCING  ON  S  -  NONE 


OBSERVATION! 


PARAMETER  PARAMETER  PARAMETER 

NUMBER  TYPE  ORDER 


mtttinmnntimtittmmutimnt  tnmn 


1  MEAN 


MOVING  AVERAGE  1 


.50644E+04 

-.42823E+00 


tttiiiilitititxilmiiiimiimiiii 


NUMBER  OF  TIME  ORIGINS  FOR.  FORECASTS  = 
NUMBER  OF  FORECASTS  AT  EACH  TIME  ORIGIN  - 
FORECAST  TIME  ORIGINS  ARE  T  = 


immtiliimimmmmmtmimim 


SDT  Tonnage  Shipped  by  MAC  to  PACAF  FORECASTS  BASE  PERIOD  40  VITH  ;5  PER  CENT  CONFIDENCE  L: 


PERIODS  AHEAD  DO.  C CNF.  LIMIT 


.  5129552E+04 
.5222385E-04 
.5222335E+04 
.  5222385E+04 
.  52223S  5S+04 
.5222335E+04 


FORECAST 


.  59091  I1E*04 
.6054430S»94 
.  5C6440QE+04 
.50544001^04 
.596440DE+04 
.  5  0  6  4  4  0  0  E + 0  4 


UP.  CONE.  LIMIT 


.66886T1Ej-04 
« 6906415E+34 
.6906415E+04 
.5906415E+G4 
.  690641  SE’-O^ 
.  5906415E*04 
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OTHER  INFORMATION  AND  RESULTS 

u  t  t  1 1 1  t  t  n  i  t  1 1  t  t  t  t  t  t  n  i  I  t  1  t  n  i  i  i  i  i  i  i  i  i  i  t  i  i  i  i  t  t  i  i  t  i  i  1  i  i  i  :  t  i  t  n i t it i n I u  n i i n i n i  i  i t  t  t i I n i i t it  n i n i i 


SAL  SUM  OF  SQUARES 


NUMBER 


RESIDUALS 


SAL  MEAN 


4 
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AUTOCORRELATION  FUNCTION 
DATA  -  THE  ESTIMATED  RESIDUALS  - 
DIFFERENCING  -  ORIGINAL  SERIES  IS  YC 
DIFFERENCES  5SLO¥  ARE 


Su;  Tarmacs 


R  DATA. 

OF  ORDER  1 


C  *  ' 

“'■•Fr 


ORIGINAL  SERIES 

MEAN  OF  THE  SERIES  =  .riCOE+03 
ST.  DEV.  OF  SERIES  =  . 44559E-04 
NUMBER  OF  OBSERVATIONS  =  36 

-.05  -.01 

.16  ,16 

-.04  .15 

.1?  .15 


MEAN  DIVIDED  BY  ST.  ERROR 


33323E-0 


SHOULD 


-"HER  THIS  SE? 
COMPARED  WITH 


15 


NOISE,  THE  VALUE 
VARIABLE  WITH  II  D 
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THE  ESTIMATED  RESIDUALS  -  SDT  Tonnaas  Shitted  tv  MSC  c:  PACAF 
CUMULATIVE  rsRIODCGRAM  .1  PROBABILITY  LIMITS  " 
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THE  ESTIMATED  RESIDUALS 
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TIME  SERIES  FORECASTING  FOR  SOT  T:nna ge  Skipped  by  MSC  to  PACAF 


DATA  -  2  -  SDT  Tonnage  Shipped  by  MSC  to  PACAF  10  OBSERVATION'S 

DIFFERENCING  ON  2  -  1  OF  GRDER  I 

tiuittniutunttutttttmttttitit  ttttutm  titniutmtti  tiimtutiinttnii 

PARAMETER  PARAMETER 
NUMBER  TYPE 

ttttmtttttitttttittttttttttt 

I  AUTOREGRESSIVE  1 

l  MOVING  AVERAGE  1 

3  MOVING  AVERAGE  1 

uttimttutttiiiuttmtnituttJiMiimtmiititmtituuntttnttti  titttttnmm  tuttt  unit  uni 

NUMBER  OF  TIME  ORIGINS  FOR  FORECASTS  =  1 

NUMBER  OF  FORECASTS  AT  EACH  TIME  ORIGIN  =  6 

FORECAST  TIME  ORIGINS  ARE  T  =  43 


PARAMETER 

ORDER 


ESTIMATED 

VALUE 


mimintit  ttnmm. 


•.26433E+00 

.74567E-01 

-.18536E-00 


SDT  Tonnage  Skipped  by  MSC  to  PACAF  FORECASTS  3ASE  PERIOD  40  WITH  95  PER  CENT  CONFIDENCE  II 


PERIODS  AHEAD 

LO.  CONF.  LIMIT 

FORECAST 

UP.  CONF .  LIMI 

1 

.28017932+05 

.37004422*05 

. 45990912*05 

L 

. 26222462*05 

. 37330332*05 

.  4853 5302*05 

3 

.22796475+35 

.  36097  61E+C5 

. 49393752+05 

4 

.2217659E+05 

.36939422+05 

. 51702252+05 

5 

. 19731045+05 

.36717162+05 

. 52653292+05 

6 

.  1812002E+05 

.  36775342+05 

. 554216~2*u5 
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This  research  was  conducted  to  analyze  the  Air  Force  Logistics 
Command  (AFLC)  forecasting  method  for  predicting  Second  Destination 
Transportation  (SDT)  tonnage  with  the  regional  flying  hour  program. 
This  thesis  had  two  main  objectives:  (1)  validate  the  current 
forecasting  method  used  for  computing  tonnage  estimates  to  derive 
SDT  budget  requests,  and  (2)  if  the  current  method's  validity  was 
not  supported,  develop  a  new  forecasting  model,  using  the  same  input 
data,  that  would  produce  more  accurate  and  reliable  tonnage  estimates. 

Analyzing  graphs  of  the  four  different  data  sets  researched  in 
this  thesis  and  then  conducting  a  statistical  test  on  the  flying 
hour  parameter  for  each  set,  it  was  determined  that  the  current 
method  employed  to  forecast  SDT  tonnage  was  statistically  invalid 
for  two  of  the  four  sets.  This  determination  was  made  due  to  the 
fact  that  the  flying  hour  parameter  changed  during  the  iterative 
regression  process  used.  This  change  implied  that  SDT  tonnage  and 
flying  hours  were  not  linearly  related. 

Box-Jenkins  (BJ)  time  series  forecasting  models  for  each  data 
set  were  built  and  provided  accurate  and  valid  forecasts.  For  the 
MAC  SDT  time  series,  the  BJ  models  were  more  accurate  than  the 
current  method.  The  BJ  models  for  the  MSC  SDT  time  series,  although 
marginally  less  accurate  than  the  current  method,  were  valid, 
whereas  the  current  method  for  these  two  series  was  statistically 
invalid. 

This  research  emphasized  the  need  for  an  accurate  model  and  an 
increase  in  the  size  of  the  data  base  used  for  forecasting.  It  was 
also  noted  that  forecasting  SDT  tonnage  requires  continual  analysis 
and  updating  to  ensure  the  model  being  used  is  appropriate  for  the 
data  being  forecasted.  Finally,  in  statistically  invalidating  the 
current  model,  this  research  has  caused  an  immediate  need  for  an 
accurate  and  valid  model.  Further  research,  particularly  with 
econometric  models,  would  prove  beneficial. 
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